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SUMMARY 


A fast, reliable computer code is described for calculating the flow 
field about a cascade of arbitrary two-dimensional airfoils. The method 
approximates the three-dimensional flow in a turbomachinery blade row by 
correcting for stream-tube convergence and radius change in the throughflow 
direction. A fully conservative solution of the full potential equation is 
combined with the finite volume technique on a body-fitted periodic mesh, 
with an artificial density imposed in the transonic region to insure 
stability and the capture of shock waves. The instructions required to set 
up and use the code are included in this report. The name of the code is 
QSONIC. A numerical example is also given to illustrate the output of the 
program. 


INTRODUCTION 

In the area of internal flows the need to design lighter, more efficient 
turbomachine components has led to the consideration of machines with fewer 
stages, each with blades capable of higher speeds and higher loading. As 
speeds increase, the numerical problems inherent in the transonic regime 
have to be resolved. These include the calculation of imbedded shock 
discontinuities and the dual nature of the governing equations, which are 
elliptic in the subcritical flow regions but become hyperbolic for 
supersonic zones. 

The present analysis combines some of the most promising transonic 
analysis techniques to calculate the flow field surrounding the cascade of 
airfoils. The full potential equation in conservative form is discretized 
at each point on a body-fitted period mesh, and a mass balance is calculated 
through the finite volume surrounding the point. Each local volume is 
corrected in the third dimension for any change in stream-tube thickness 
along the stream tube. The nonlinear equations for all volumes are of mixed 
type (elliptic or hyperbolic) depending on the local Mach number. The final 
result is a block-tridiagonal matrix formulation involving potential 
corrections at each grid point as the unknowns. The residual of each 
equation is equal to the net mass flux through a volume. The resulting 
system of equations is solved along each grid line by an SLOR procedure. 

At points where the Mach number exceeds unity, the density at the forward 
(sweeping) edge of the volume is replaced by an artificial density (ref. 1). 

The capabilities of this FORTRAN program, beyond those of previous 
cascade analysis methods (ref. 2), include the ability to calculate through 
weak shocks (peak relative Mach number less than 1.4) and completely around 
both leading- and trai 1 ing-edge regions of a blade. The blade shapes in the 
leading- and trai 1 ing-edge regions are not restricted to circular arcs. 

The flow field to be calculated can be axial, radial, or mixed and can 
be rotating or stationary. The geometry of the stream surface can be a 
two-dimensional planar cascade or axi symmetric with varying channel 
thickness and radial position. There are, however, several situations that 
cause unsatisfactory results because of the inviscid flow assumptions. 

Sharp leading edges at high incidence and high-Mach-number turbine blade 
trailing edges with substantial deviation will both cause large velocity 
peaks on the blade as the inviscid flow tries to negotiate the large angle 
from stagnation point to suction surface. In addition, if the passage is 
nearly choked, the code may have difficulty converging. 


More details on the methods used by the code for grid generation are 
discussed in reference 3. In the following sections an overview of the 
program's capabilities is given, and an outline of the solution procedure is 
presented. A list of the necessary input, an interpretation of the output, 
and the basic flow logic of the code are given. Sufficient information is 
presented to allow a user to make use of the program. A detailed test case 
is also included for ease of verification. 


PROGRAM DEFINITION 

The QSONIC code was developed to obtain the transonic potential flow 
distribution on a blade-to-blade stream surface of a stationary or rotating 
turbomachine blade row. 


Capabi 1 ities 

The following assumptions were made in developing the analysis and the 
code: 

(1) The flow is inviscid and adiabatic. 

(2) The flow relative to the blade is steady. 

(3) The fluid is a perfect gas with constant specific heat. 

(4) The flow is isentropic and any discontinuities (shocks) are weak 
enough to be approximated as isentropic jumps. 

(5) There is no velocity component normal to the stream surface 

(fig. 1). Flow is restricted to the two components within the stream 
surface: the streamwise direction (vector combination of axial and radial), 

and the tangential direction. 

(6) The flow relative to a frame fixed in space (absolute velocity) is 
completely irrotational, and the flow relative to the blade can be 
represented as the vector addition of this absolute flow field minus the 
wheel speed in the tangential direction. 

Assumption (6) is sufficient to allow the absolute velocity to be 
represented as a gradient of a potential field. 

The flow field to be calculated can be axial, radial, or mixed. The 
compressor or turbine blade row can be fixed or rotating. As a result of 
assumption (4) the peak local relative Mach number should be 1.4 or less. 

The geometry of the stream surface can be a two-dimensional planar cascade 
or an axisymmetric configuration whose stream-sheet thickness and radius 
vary with distance along the stream surface. The extent of this surface 
that corresponds to one blade passage, with the blade in the center (see 
upper portion of fig. 1), is mapped into a rectangle for numerical 
simplicity by one of the two mesh generators (ref. 3) available in QSONIC. 

Figure 2 gives major features of the solution process, including grid 
generation and refinement. There are essentially two modes in which the 
program can be operated: (1) One grid generation and a flow solution can 

be combined in a single execution. (2) A separate run is first made to 
generate and store a grid, and then the flow solver is used sequentially on 
several different meshes that are based on the stored grid (mesh 
refinement). Mesh refinement normally gives about a 30 percent gain in 
convergence speed over what can be attained by using the finest (stored) 
grid exclusively to meet a given convergence criterion. 
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Figure I, - Schematic of stream-channel surface. 
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Figure 2. - OSONIC solution process. 

















Limitations 


Flow Solution Method 

Although the program will perform well for a wide variety of turbine and 
compressor blading and flow conditions, there are several situations that 
cause unsatisfactory results because of the potential flow assumptions. 

Since the flow is to be inviscid and irrotational , sharp leading edges at 
high incidence (more than a few degrees) will cause large velocity peaks on 
the blade surface as the inviscid flow tries to negotiate the large angle 
turn from stagnation point to suction surface. Such peaks should be 
discounted. The remainder of the flow field will still be reliable since 
these peaks seldom poison more than a few grid points in their vicinity. 

This situation arises most often with thin-blade compressor geometries, 
where the adverse flow gradient on the suction surfaces tends to discourage 
the velocity fluctuations from propagating very far or reflecting off the 
downstream boundaries. 

The conditions are not so favorable at the trailing edge of 
high-Mach-number turbine blades, where a more serious problem can arise. 

Here the flow is at its highest velocity, the blade has its minimum 
thickness, proximity to neighboring blades is enhanced (in the area of the 
passage throat), and the flow gradient is favorable for carrying errors and 
velocity peaks toward the downstream boundary and further reflection. 
Velocity spikes that arise here can be very severe, propagating across the 
narrow passage and eventually choking or poisoning the entire region. The 
safest method for using the code on geometries such as this is to approach 
the actual conditions slowly. Begin with a coarse grid (no more than 14 
cells across the passage) and slowly increase the upstream Mach number or 
vary the downstream flow angle (away from minimum deviation) in the 
neighborhood of the desired values. If more than three points on the blade 
surface exceed Mach 2, the solution is not reliable. 

A third and final area of difficulty for the code is choking flow since 
the mass flow is determined by the upstream boundary conditions and global 
continuity is assumed. But whether a particular shape of passage is choked 
or not can only be determined by executing the code. The program may run to 
completion for nearly choked passages (maximum mass flow), but the solution 
would be unreliable if there were an area of Mach 1.1 extending more than 
halfway across the passage. 


Grid Generation Method 

The generator referred to as the interpolation scheme (ref. 3) is best 
applied to compressor or turbine blading of low to moderate camber and 
thickness. ^This mesh generation scheme has no problem with large stagger 
angles (=75*) nor with high solidities. Cases with solidities of 1.45 have 
been attempted with good results. However, if the turning angle or 
thickness of the blade is large, while at the same time stagger and solidity 
are near the limits mentioned above (i.e., stagger plus turning =124”), the 
grid lines may become crowded on the suction surface. This is common for 
turbine blade geometries. The interpolation scheme does not perform well 
for blades with very small leading- or trai 1 ing-edge radii (knife edges). 

It does have the advantages of allowing some control over placement of grid 
lines by concentrating them near the blade surface or in the trai ling-edge 
region and also of blade coordinate smoothing. 



(a) Interpolation scheme grid. 



(b) Electrostatic analog grid. 

Figure 3. - Coarse mesh for turbine cascade. 


The electrostatic analog scheme (ref. 3) was designed to handle thick 
turbine blades of high camber, thickness, and solidity, so it is general 
enough to handle any possible blade geometry, including knife-edge 
profiles. The disadvantages are the lack of control over grid line 
concentration and no smoothing. Figure 3 shows results of the two 
generators for the same turbine cascade. Both grids shown were successful 
in subsequent flow solution runs. 


Program Operations 

The FORTRAN code for QSONIC was developed on an IBM 3033/370. The 
accuracy for the present version of QSONIC is limited to single precision 
(for this machine, seven decimal digits). The program requires 58 000 words 
or core storage to execute a separate geometry generation run with the 
interpolation scheme. This rises to 63 000 words for an electrostatic 
analog grid generation. An execution that combines either of the grid 
generators with a flow solution requires 210 000 words, provided that the 
common blocks (ENTIRE, 6E0M, and GE0M2) containing the mesh coordinates and 
parameters are loaded separately and the grid generator codes are overlayed 
with the flow solver module. If the entire code is loaded simultaneously, 
the maximum core size needed is 240 000 words. If common block HIBALL is 
eliminated (used for metric storage to avoid recalculation), this reduces to 
90 000 words. In all applications the BLOCK DATA module should be the first 
one loaded. 
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SOLUTION METHOD 


Governing Equations 

The flow field relative to a frame fixed in space (the absolute system) 
is assumed to be irrotational , isentropic, and inviscid. Thus the 
continuity equation 


If + V • (pV) = 0 


( 1 ) 


where V is the absolute velocity, can be rewritten with V being the 
gradient of a potential field 


Y = Vip 


and 






(pVtp) = 0 


J 


( 2 ) 


But this field is unsteady because of the passing blade. If we assume that 
the coordinate system is fixed to the blade, continuity takes the desired 
steady form 


V • (pW) = 0 


where 


TI = V - ft X r 


j 


(3) 


which is just the absolute velocity minus the wheel speed in the tangential 
direction. 

It is required that equation (3) hold at each point in the flow field. 
As suggested by several others (refs. 4 and 5), equation (3) is integrated 
over the element of volume at each point. Making use of the divergence 
theorem gives 



W dAg = 0 


(4) 


where n is the outward normal to the area element dA, Y the element 
volume, and Ag the surface area of the element. 

Approximating the integral as the sum over the four element faces shown 
in figure 4 gives 
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If we differentiate equation (5a) to form 6R*^, equation (6a) becomes 


6^ 


■Z( 




^ K . 
n6p aAp 


p n 


6 ^ 


aA, 




p W 


n aAc 


But for isentropic flow 


2 - 




and for this potential flow 


(7) 


(8) 


W = V(p - h X r (9) 

Then can be expressed as a function of the local geometry and the 

potential values for <p. From equation (8), since p is a function of W, 

6p*^=|^6w'^ (10) 

From equation (9),^after transformation to the computational (n,5) space, 
for 6W = 6ui + 6VJ, 

= 3(<Scp)*^ 3T1 ^ 3(6ip)*^ 3C 

3n 3X 35 3X 

and 

_ 3('S9)^ 3n ^ 3(6(p)*^ 3C 

3n 3y 35 3y 

where 

, K , K+1 Kv 1 

6(p = (ep _ q, ) _ 

0) 


( 11 ) 


( 12 ) 


(13) 


is the unknown correction to the local potential field, and u is an 
overrelaxation factor to accelerate convergence. 

If we substitute equations (8), (10), (11), and (12) into equation (7) 
and recall that the area aAe; of each face of the element is just the 
product of its length i and stream-channel thickness h, we obtain 
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/ K . K \ 1 / j/ . K { 
y '' ®yj 3W K P ®x V 


a(6<p)*^ In + a(gy)*^ il\ 
an ax ac ax / 


fu'^e + v'^e ^ ^ + p'^e iH + il) i &h 

\ X '' ®y/ aW yK ^ y \ an ay ac ay/ [ 


=-s 


p*^(ue^ + ve^)'^Ah (14) 


where the normal vector to each face is 

n = e^T + e^j (15) 

The right side of equation (14), Ri^j, is the residual net mass flow 
out of the finite volume surrounding the grid point (i,j). Each volume 
element is defined by the four surfaces parallel to and located halfway 
between neighboring grid lines, as shown in figure 4. Several schematic 
views of the quasi-three-dimensional stream surface are presented in 
figure 1. The thickness of the element (i.e., relative local stream-tube 
height) and its radial position are obtained from specified distributions. 
Height is set to unity for planar two-dimensional flow. Variations in the 
radius r of the stream tube are accounted for by scaling tangential re 
coordinates y with a ratio of the local radius r to a given upstream 
radius r[_E. Then over the area Ag defined by the four surfaces of the 
volume element. 




K ttK , 

i±l/2,j‘^i±l/2,jS-±l/2,j^i±l/2,j 


1^ \ 

Pi,j±l/2^i,j±l/2S,j±l/2'^i,j±l/2j 

where h is the_local stream-tube height normalized to unity far upstream. 
The velocities W are formal to the element boundaries and located at the 
boundaries. Velocity W and density p are calculated from values of 
potential at the previous iteration. 

The summation shown on the left side of equation (14) proceeds in the 
same manner as that for R. Quantities are calculated at the center of each 
face (fig. 4), (i, j + 1/2), (i, j - 1/2), (i + 1/2, j), (i - 1/2, j), and 
finite differences are formed by using potential correction values at 
neighboring points. As an example. 
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an 


, K 


- 6<p 


111 


i,j+l/2 


n,- 


j+1 


- n. 


(17) 


and for cross terms 


3(6«p) 


ac 


i,j+l/2 


= 1/2 


+ '■'in 


i+1 


- t 


1 


K-1 , K \ 

J 


(18) 


If we make use of the relations and similar expressions for the other 
faces and isolate all unknown potential corrections on the left side, 
equation (14) can be expressed in the form 


K 


K 


^i,j+l'^‘‘’i,j+l ^i,j®‘‘’i,j 


+ A, 


. . , 69 '^ . , = . 

i,j-l ’i,j-l 


K 


(19) 


where 






= -^i 


+ G 


/ K K , K , K-l\ 

(9 ,P + I 


(20) 


By writing equation (19) for each point on a grid line of constant c 

(j = 1,J), a tridiagonal matrix system is formed that is then solved by the 

standard Thomas algorithm (ref. 6) for the potential corrections 69 *^. 

New values of the field potential are obtained by using equation (13) after 
all lines are solved (i = 1,1). 

The coefficients in equation (19) are 
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Sweep Strategy 


The grid lines are solved in order, beginning with the line that goes to 
downstream infinity and then proceeding clockwise around the blade shape. 
This order results in a solution that sweeps against the flow direction on 
the bottom surface of a blade. It became evident, when symmetric supersonic 
flows were successfully solved, that direction of sweep relative to the flow 
is not critical. The apparent reason is that additional damping is applied 
by the elliptic-like nature of the artificial density when it is used in 
combination with the sweep strategy. More details on the effect are given 
in the next section. 


Artificial Density 

In equation (7) there is no explicit artificial viscosity term necessary 
for stability in supercritical flow regions and for smearing shock jumps. 

As suggested by Hafez et al. (ref. 1), an implicit viscosity term is 
incorporated in the density. This modified density is used instead of the 
physical density in equation (7). Subsequent numerical experiments 
indicated that modifying the densities on the left side of equation (7) at 
each iteration level was not necessary. Thus, to speed the calculation, 
there is a parameter (IT) that controls how often to update the densities. 
The form of the artificial viscosity modification is (ref. 1) 


p = 


p - 


lS 


( 22 ) 


where y is the switching function to eliminate the density modification in 
subsonic regions (M < 1) 



(23) 


where M is the Mach number relative to the moving blade row and ps aS 
is the total differential component of the density in the flow direction S 
defined by the contravari ant velocity components a and b. Then 


MPS = Pg - p^ An + - p^ AC 


(24) 


where 





i£ 

9C 
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The switching function was split into two terms so that the effect of each 
component of the density change could be even more efficiently gaged in 
areas of high flow gradient. In addition, an upwind-averaged Mach number 
was used: 



where 


M 


av. 


= T (M. 


i+l/2.j 


'^i+l/2,j±l^ 


Wjj = max 




( 25 ) 


where 




i+l/2.j 


+ M 


i+l/2±l. 




J 


Similarly, upwind differencing is used for approximating and 

The physical densities used to calculate p^ and p^ are 
updated as each line is reached. Implicit in this arrangement is the 
creation of a timelike cross derivative of density pg^ when sweeping 
against the flow direction. This cross derivative is created because the 
density at the current line has been calculated by using the values of 
potential and velocity at the latest iteration (i.e., artificial time step), 
but the values at the next upstream line, where the upwind differences are 
formed, remain at the previous time step until they in turn are solved. 
However, when sweeping with the flow, the upstream line from which the 
density gradients are formed has already been solved and updated. Thus this 
additional damping due to the time lag of densities is not present when 
solving with the flow, and the artificial compressibility reduces to the 
original space lag only. 


Boundary Conditions 


The boundary conditions applied to the solution system include a 
specified upstream axial mass flow, in the form of an oncoming Mach number, 
plus incidence and stagger angles. Using global continuity. 


^out^out^out ^in*^in^in 


( 26 ) 
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where A is the flow area normal to the velocity u (fig. 1). 

The specified flow angle at the downstream boundary, together with a 
simple iteration procedure on the isentropic density ratio, gives the 
downstream Mach number boundary condition. For a stagger angle x and 
upstream relative angle of attack equation (26) can be expressed as 


Y+1 


out 


1 + 


Y — 1 m2 


A. cos a. 

F<«out> 


out> 


'out 


‘out 





where a = 6 + x, n is the rotational speed, and r is the radius of the 
stream channel. Using a straightforward Newton iteration on equation (27) 
yields Mq^^. From this the relative speed ratio is available as 



Through trigonometric manipulation of the upstream and downstream velocity 
triangles, by using the input values of stagger, relative incidence, and 
rotational speed, the absolute flow angles are obtained. The values 
of the velocities at volume elements that touch the far-field boundary are 
then calculated by using equation(26) and assuming constant angular momentum 
in the far field. Thus 


puA 


‘^in^in^in 


(29) 


with 


rv = -^ = Constant (30) 

3 w 

Since density is a function of velocity, equation (29) can be rearranged 
for a Newton iteration procedure to obtain u at each border element point. 
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f(u) = u - 


0 


(31) 


p(u,v)A 

Periodicity is maintained in the tangential (y) direction by enforcing the 
relation 


9 (x, y + s) = ip(x,y) + C 


(32) 


where s is the blade spacing and 
,y+s) 




= Vin^in 


(33) 


On the grid line defining the body, a standard zero-normal flow condition is 
imposed. 


(vq> - orj) *0 = 0 


(34) 


DESCRIPTION OF SUBROUTINES 

Figure 5 illustrates the calling relationships between the main module 

and its subroutines. 

Executive 

MAIN QSONIC Reads case control namelist for each job. Calls 

geometry setup, stream-channel data supplier, and flow 
solver for each of several grids requested (mesh 
refinement). References final output generator after 
last mesh is solved. 

BLOCK DATA Initializes all data values for variables that are in 

common blocks. Nonexecutable. 


Geometry Generation 


WRAPUP 


PJGRID 


Sets up geometry. For a new blade, reads bulk 
namelist of body coordinates, normalizes all 
dimensions to aerodynamic chord, and calls the mesh 
generator that was requested. For geometries already 
developed, reads file containing grid of maximum 
fineness and selects grid lines at proper intervals to 
create the coarser grid requested. Rearranges data 
arrays so that (1) the first row of each matrix 
contains the line going to downstream infinity, 

(2) the aerodynamic chord lies along the x axis, and 

(3) the origin is at midchord. 

Generates a mesh by using the electrostatic analog 
method (ref. 3). 
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MAIN 
OSONIC 
5:R; 23: B; 
13: B 


RBIN 

2:R 


FILRB 


GRIDDL 


WRAPUP 
5:R; 23: R; 
18:B,W; 13:W 


WHEEL 


FILPHI 

18:B,R; 8:B,W; 
9:R,B 


PJGRID 


uso 

NIC 

KEEPER 
7:B; 13:B,W 



SPEEDS 


ALTERP 


MMDELK 


BADZAC H 


OWOW H 


TOSTAG 


FUDS 


u 


JELF 


SURFUP 


O' 


FZTRP 

h 

1 

XYCALC 


SGEN 


T 


DZTRP 

I5I 

DNTRPC 


SPGEN 


->-[ FNTRP 




FNTRPC 




T 

XNTRP 

IZTRP 

RZTRP 


CHLSKY 


NRRTS3 


SIZE 


ORDER 


BNTRP 


'5' I/O unit referenced 
B Rewind 
R Read 
W Write 


FNTRPA 


TLU 


LIMIT 


Figure 5. - Calling relationship between modules. 
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CHLSKY 

GRIDDL 

SURFUP 

T0STA6 

ALTERP 

FUDS 

OWOW 

BADZAC 

JELF 

MMDELK 

(lOPT.ARG 

WHEEL 

Interpol a 
XYCALC 


Uses Cholesky's method to solve a linear system of 
equations (ref. 7). 

Generates a mesh by using the interpolation scheme: 

(1) Performs analytical transformation steps 

(2) Calls module to incorporate blade shape 
within flat-plate grid 

Stretches the coordinates of the flat-plate mesh so 
that the desired blade surface becomes the innermost 
closed curve (grid boundary). 

Transforms the flat-plate cascade at zero stagger into 
the circle plane. Rotates the circle to the desired 
stagger and inverts the transformation back to the 
physical plane. 

Interpolates on tabular values of solidity and stagger 
for the location of the positive and negative infinity 
points in the circle plane. 

Calculates one of the functions needed (defined as T 
in eq. (1), ref. 3) for transformation of the 
flat-plate cascade into the rectangle in computational 
space. Calls routine to calculate Jacobian elliptic 
functions. 

Calculates the function defined as n in equation 
(1), reference 3, needed for the transformation to a 
rectangle. 

Calculates the function defined as r in equation 
(1), reference 3. 

Calculations Jacobian elliptic functions SN, CN, and 
DN. 

Computes double-precision value of complete elliptic 
lER) integral of the first kind, for d<p/SQRT. 

(1. - ARG • SIn2<p). IOPT = 1, 0 < ARG < 1. 

Rotates x,y coordinates to the stagger angle. 


ion and Smoothing 

Smooths the curve defined with the input points by 
generating additional points along the curve. The 
input is used to develop double three-point 
interpolating polynomials in successive regions along 
curve. Polynomials are then used to suggest points 
that can be tested against the spacing criteria as 
defined in the SPGEN listing. 
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SGEN 

SPGEN 

LIMIT (I,J,K) 

DZTRP 

FNTRP 

FNTRPA 

DNTRPC 

FNTRPC 

TLU (TABLE, 

ARG,N,I) 

FZTRP (S,F,A,N) 

XNTRP(S,F,FO,N, 
ANS), with entry 
points IZTRP(S,F, 
FO,N,ANS) and 
RZTRP(S,F,FO,N,ANS) 

BNTRP 


NRRTS3 

SIZE (X,I,N), with 
entry point ORDER 
(X,I,N) 

Stream-Channel Data 
RBIN 


Computes line integral of the distance along the curve 
defined by the array of input point pairs. 

Generates a table of points along the line integral as 
widely spaced as possible (see SPGEN listing for 
details) . 

Limits J between I and K, i.e., if J < I, J = I; if 
J > K, J = K. 

Evaluates complex derivative for three-point 
interpolation. 

Evaluates real function for three-point interpolation. 

Evaluates first coefficients for three-point 
interpolation. 

Evaluates final coefficients of derivatives for 
three-point interpolation. 

Evaluates final function coefficients for three-point 
interpolation. 

Table lookup of I, such that TABLE(I) < ARG < 

TABLE (I + 1) 

Interpolates on line integral array S, corresponding 
to point pairs in F, to find A on S. Returns matching 
point of F. 

Inverse of function FZTRP, which interpolates on 
complex array F to find the imaginary part (if 
IZTRP is referenced) or the real part (if RZTRP is 
referenced) that equals FO. Returns the matching 
value from S as ANS. 

Calculates polynomial coefficients needed for XNTRP 
interpolation. 

Finds real roots of a cubic equation. 

Loads the array I with indices from X array in the 
order of increasing X value. Rearranges X according 
to pointer array I. 


^reparation 


Reads in the stream-channel position and thickness 
arrays from a file with the identical format as the 
input of the TSONIC code (ref. 2). 
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FILLRB 


Fills the arrays of radius and stream-tube thickness 
that correspond to the points in the solution grid by 
interpolating on the data from RBIN or by assuming a 
linear variation in the x direction. 


Flow Solution 

USONIC Calculates the downstream boundary conditions, 

initializes the potential function, and relaxes the 
compressible potential flow equation on a given mesh 
until tolerance is reached. Calculates flow 
velocities on. the blade surface when tolerance is 
reached. 

FILPHI Based on the coarse mesh solution for the potential, 

interpolates values for the next, finer mesh. 


Output Generation 

KEEPER (an entry Calculates and prints velocities, densities, pressures, 

point in USONIC) and Mach numbers at all points in mesh. Records all 

results on file for later graphic display. 


OPERATING INSTRUCTIONS 

The information required to operate the mesh generators and flow solver 
consists first of a two-dimensional description of the blade shape in the 
stream surface. This is in the form of pairs of (x,y) points on the surface 
(fig. 6), together with parameters that describe the cascade layout (chord, 
stagger angle, and pitch). A group of parameters describing the desired 
density of the mesh lines completes the input for a geometry generation. To 
calculate the flow field, the upstream and downstream flow conditions, the 
rotational speed, the convergence criteria, and a schedule of meshes to be 
used should be input. If quasi-three-dimensional effects are to be 
included, a file containing a description of the stream channel's radial 
thickness and position as a function of distance along the stream surface is 
needed. The format for this file is identical to the input form for the 
code described in reference 2. 

Output of QSONIC consists of printed listings, a formatted file with all 
generated mesh points, and a formatted file with final velocity components, 
Mach numbers, and pressures at all mesh points for use in graphic displays. 
The printed listing contains repeated input data, generated mesh coordinates 
on the blade, periodic progress reports on the flow convergence, and a list 
of the final velocities, pressures, and densities on the blade surface for 
each grid included in the schedule of solution meshes. 
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(b) In meridional surface. 
Figure 6. - Geometric input variables. 
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Input/Output Files 


I/O unit 


Usage 


2 


File containing stream-channel data in TSONIC input 
format 


5 Standard card input; NAMELISTS PARAMS, INSTUF 

6 Standard printed output 

8 Restart save file (read as unit 9 when restarting) 

9 Restart input file, containing last successful print 

of potential values 


13 For mesh generation runs, coordinates of all mesh 

points are written here for future solutions (415, 
followed by all 6E13.5). For flow solutions, x, y, 
velocities, pressures, and Mach numbers are recorded 
for later graphic display (415, 6E13.5). 

18 Temporary storage 

23 Previously developed mesh coordinates are read in from 

this file (= 13, from mesh generation runs). 


Input 

All input for any run will fall into the following categories: 

(1) Logical and case control parameters (NAMELIST PARAMS, on unit 5) 

(2) Bulk data input: 

(a) For geometry generation runs (NAMELIST INSTUF, on unit 5) 

(b) For flow solution runs: 

(i) Mesh point storage file (unit 23, equals output from geometry run) 

(ii) Stream-channel data file (unit 2) 

(iii) Restart input file (unit 9) 


Logical and Case Control Parameters 

Except for TITLE, the format for all of these variables is in NAMELIST 
form. The NAMELIST is PARAMS. If a default value is listed as none, a 
value of that parameter must be input. N/\MELIST variables can be entered in 
random order. 

Parameters applying to all runs. - The following parameters apply 
to all runs: 
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Name 


Type 


Default 


Description 


TITLE 

NOFLOW 

RESTAR 

REMESH 

MS 

NOZES 


MESHl 


Alphameric, 

(format, 

20A4) 

NONE 

A one-line name for the case being run; 
TITLE must appear on the first line of 
unit 5. 

Logical 

FALSE 

=TRUE if run is to stop after generating 
a mesh (a geometry generation run). 

Logical 

FALSE 

=TRUE if run is a continuation of a 
previous effort. 

Integer 

4 

Index in MS and NOZES arrays of the 
first grid to be used when restarting. 
Only flow solution runs may be restarted. 

Integer 

NONE 

Array (dimension 10) of values (maximum 
value, 25) for the number of grid lines 
that will enclose the blade. Twenty-five 
is a satisfactory number for blades of 
solidity near unity. As solidity 
increases, the maximum in MS should 
decrease to avoid crowding. Twenty-five 
allows five extra lines for high-camber 
blades (see CAMPER). The absolute limit 
is 30. 

Integer 

NONE 

Array (dimension 10) of values (maximum 


value, 49) for the number of grid lines 
radiating from the blade on one surface, 
starting with the line going to upstream 
infinity and ending with the line going 
to downstream infinity. Forty-nine 
radiating lines and 25 enclosing lines 
gave fairly square grid elements for 
solidities near 1. Maintaining similar 
aspect ratios between successive grids 
during a solution is a conservative rule 
to avoid convergence problems. If a 
value of NOZES is greater than zero, 
a mesh with that many lines will be 
developed and stored. If the value is 
negative, it is assumed that this mesh 
already exists on file and will be read 
in. 

Integer 1 Index in MS and NOZES of first grid to 

be calculated (geometry generation) 
and/or used for flow solution. MS and 
NOZES offer a choice of grids. For 
geometry generation runs, MESHl selects 
one of these to be developed and stored, 
provided that NOZES (MESHl) > 0. 
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MESHN Integer NONE Index in MS and NOZES of last grid to be 

calculated and/or used for flow solution. 
Normally, for geometry generation, 

MESHN = MESHl. Note: The purpose of a 

geometry generation run is to develop a 
grid of the maximum fineness desired. 
Subsequent flow solution runs can then 
solve the case on all grids listed in MS 
and NOZES between index MESHl and MESHN 
by selecting an increasing number of the 
grid lines available on file ("grid 
refinement"). Grid refinement normally 
gives about a 30 percent gain in 
convergence speed over what can be 
attained by using the finest (stored) 
grid exclusively to meet a given 
convergence criterion. A schedule of 
three increasingly dense grids is 
usually sufficient. For flow solutions, 
MESHl normally points to the coarsest 
mesh. To insure equal spacing of grid 
lines over the whole mesh, it is often 
helpful to require that 

NOZES(I) - 1 
^(I-MESHI) 

be an integer, where MESHl £ I £ MESHN, 
and K = 2, 3, 4, 5, etc. (i.e., K = 2 if 
number of grid lines is doubling between 
successive grids). A similar result 
holds for the values of MS. The logic 
of the code is much more general than 
these restrictions would indicate. 
However, whether the resulting grids 





would be suitable for flow solution is 
not certain. 

LAMDAO 

Real 

NONE 

Stagger or setting angle of the blade 
row, in degrees measured from the 
throughflow (x) direction to the blade 
chord line (clockwise negative). See 
figure 6. 

CHORD 

Real 

NONE 

True chord of the blade, in same units 
as blade coordinates. 

S 

Real 

NONE 

Blade spacing in y(re) direction, 
same units as coordinates. 


Parameters required only for geometry generation. - There are two grid 
generators available, as decribed in reference 3. The electrostatic analog 
is applicable to any blade shape and to any value of stagger or turning. 

The interpolation scheme is good for most blades except knife edges (i.e.. 
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edge radius < 0.005 chord) or when the sum of stagger plus turning angles 
approaches 120°. Grid lines can be concentrated near the trailing edge or 
the periodic boundary or around the blade. No concentration is available 
when using the electrostatic analog. Refer to figure 6 for illustration of 
geometric input variables. 


Parameters required for either grid generator: 


Name Type 

NED Integer 


KN Integer 


Default Description 

NONE Total number of body definition 

coordinates to be input; one table of 
points with first point repeated as last 
point. See the section Bulk data for 
geometry runs . If NED < 0, no blade is 
read in and a flat-plate cascade is used. 

NONE Flag to indicate which grid generator is 
to be used. For electrostatic analog, 

KN = 0 is input. For interpolation 
scheme, KN = number of body points on 
upper surface (from minimum to maximum 
X points, inclusive). 


For the electrostatic analog generation, no other control parameters are 
needed. Refer to the section Bulk data for geometry runs for instructions 
on preparing the body points. Note: For this generator, the value from 

NOZES that is used must be odd. 


Parameters applying only to interpolation scheme: 


RLE 

Real 

NONE 

Leading-edge radius of blade. Same 
units as coordinates. See figure 6. 

RTE 

Real 

NONE 

Trai 1 ing-edge radius of blade. Same 
units as coordinates. 

THETL 

Real 

0 

Camber angle of leading edge, deg, 
measured from the aerodynamic chord to 
the line tangent to the mean camber 
curve at the leading edge. Clockwise is 
positive. 

THETT 

Real 

0 

Camber angle at trailing edge, deg, 
measured from chord line to the line 
tangent to mean camber curve at the 
trailing edge. Clockwise is positive. 
Specification of THETT and THETL gives 
better interpolation results for thin 
blades with high camber. 
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CAMPER 

Integer 

6 

For blades whose chordline lies outside 
the blade, extra grid lines surrounding 
the blade are needed to interpolate the 
blade position. The truncated value of 
MS( ) /CAMPER is added to MS( ). The 
maximum allowed value of MS( ) + MS( )/ 
CAMPER is 30 grid lines (see MS). 

STABAC 

Real 

0.999 

For test cases where no blade shape is 
to be input (NED < 0) and if STABAC = 
0.999, a flat-plate grid is developed. 

For 0. < STABAC < 0.999, one of the 
surrounding grid lines will be designated 
as a body shape and a mesh developed 
around it. As STABAC approaches zero, 
the grid line designated is ever closer 
to the periodic boundary. 

CHOP 

Real 

0.99 

Fraction of distance from midchord to 
either end beyond which all lines that 
intersect the flat plate near the 
leading or trailing edge will receive 
special handling to calculate the point 
where they intersect the blade shape. 
This is often necessary for blades with 
small leading- or trai 1 ing-edge radii 
(less than 2 percent chord). Normally, 
0.9 < CHOP < 1.0. 

SMOOTH 

Logical 

False 

= TRUE, for automatic addition of more 
blade definition points in the region 
of both leading and trailing edges. 
Recommended for all blades except cusps 
and wedges. 

LEONLY 

Logical 

False 

= TRUE, for smoothing about leading 
edge only, for the case where only the 
trailing edge has a cusp or wedge. 
Smooth must also be .TRUE. 

SLPl 

Real 

1. 

If SLPl < 0., concentration of grid 

SLP2 

Real 

2. 

lines is desired in some region of the 

SLP3 

Real 

1. 

mesh. These four parameters control the 

SLP4 

Real 

1. 

amount and location of the concentration. 
If SLPl < 0., its absolute value and 


that of SLP2 refer to the end-point 
slopes on a fifth-power polynomial 
curve of the new values of 5 (referred 
to as 5 ' in fig. 7) as a function of 
the original values of where 
5 -lines are those that radiate from 
the blade; but if SLP2 is also less than 
0., its absolute value is interpreted 
as ds'/dclpK* where FK' is the value 
of 5 at the trailing edge. A small 
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lfSLP2<0: 




(b) Surrounding contours. 

Figure 7. - Polynomial curves used for concentration of grid lines. 


value of |SLP2| in this case results in 
a flatter curve in the trai ling-edge 
region, thus concentrating more values of 
c' in that region. If the concentration 
is not desired at the trai-ling edge, the 
value of F (the variable name for F is 
EF), between 0. and 1., which determines 
the point of concentration on the blade, 
can be adjusted by manipulating input 
values of SLPl and SLP2 as shown in the 
following table. Exact control is not 
available in this version of the code. 
Hence concentration should be used 
cautiously. 
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I 


Near line 
going to 

Value 
of EF 

Input 

SLPl 

Input 

SLP2 

+00 

0 

.25 

-1.533 

+17. 


.5 

-1.667 

+5. 


.75 

-2.1428 

+2.777 

— oo 

.99 

-26. 

+2.02 


SLP3 and SLP4 refer to the end-point 
slopes on a similar curve of new values 
for n (n* in fig. 7) as a function of 
the old n values. SLP3 and SLP4 
would normally be less than 1., since 
concentration near the end points 
(periodic line and body) is desired 
for n. 

Parameters required only for flow solution runs . - The following 
parameters are required only for flow solution runs: 


Name 

Type 

Default 

Description 

MINF 

Real 

NONE 

Mach number at upstream boundary. For 
rotors this is the relative Mach number. 
If the downstream Mach number is known 
to be supersonic, input MINF with a 
minus sign. 

BETAl 

Real 

NONE 

Flow angle at upstream boundary, deg. 
For rotors this is the angle in the 
relative frame. BETAl is measured 
from the aerodynamic chordline to the 
direction of flow. Clockwise is 
negative. 

BETA2 

Real 

NONE 

Flow angle at downstream boundary, deg, 
measured from aerodynamic chordline 
to direction of flow. Clockwise is 
negative. For rotors this is the angle 
in the relative frame. 

GAMMA 

Real 

1.4 

Ratio of specific heats 

TOLS 

Real 

NONE 

Array (dimension 10) of tolerances 
corresponding to MS and NOZES. Each 
grid solution will proceed until its 
TOLS value is satisfied. There are 
three forms of input permitted: 
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(1) -1. < TOLS(I) < 0.: Proceed 

until the relative circulation error 

I (Ccalc “ 9s^act) /^ex^ct I 

than |T0LS(I)|, This is appropriate only 
for lifting (nonsymmetric) cases and is 
indicative only of convergence in the 
trail ing-edge region. ,T0LS values 
between -10“-^ and -10~° are typical 
for most grids with the input value 
approaching zero as the mesh becomes 
finer. 


(2) 0. < TOLS(I) < 1.: Proceed until 

the average relative change in potential 
on the mesh is less than 

|T0LS(I)|. This is a more rigorous 
criterion but slower to be met because 
of the averaging of all grid points. 
Values typically should be between 
10~3 and 10“^. There is no benefit 
to lingering on the coarser grids much 
after the average correction has leveled 
off . 


(3) 1. < TOLS(I): Proceed until the 

number of iterations equals TOLS(I). 
Regardless of the TOLS input, each grid 
solution will terminate after 300 
iterations if the TOLS criterion has not 
yet been met. 


OVEREL 

Real 

1.5 

Overrelaxation parameter, 1. < w < 2., 
for all subsonic points in the mesh 

UNDERL 

Real 

1. 

Underrelaxation parameter between 0.5 
and OVEREL. During the first few 
iterations on each grid to is changed 
gradually from UNDERL to OVEREL. 
Default is usually adequate. 

SUPREL 

Real 

1. 

Overrelaxation parameter applied to 
supersonic points. Maximum value, 2.0. 

NOWREL 

Integer 

20 

Number of iterations on each grid before 
subsonic overrelaxation is applied. 
Normally 2 < NOWREL < 20. 

NOTYET 

Integer 

2 

Number of iterations on each grid before 


compressibility effects (p/pin ^ !•) 
are allowed. SUPREL is also turned on 
after this many iterations. Normally 
2 < NOTYET £ 10. If NOWREL > NOTYET, 
UNTTERL is not used. 
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TEGARD 


Real 


DAMP 

CII 

IT 

ALLOUT 
QUASI 3 
NSTRM 

RINF 


WAKE 


Real 


Real 


Integer 


Logical 


Logical 


Integer 


Real 


Real 


2.0 Default is adequate, except for 
high-speed blades with a wedge or small 
radius in the leading- or trai 1 ing-edge 
region. Increases damping in edge 
regions. Use cautiously in combination 
with DAMP to avoid instability. 

1. < TEGARD < 10. 

1.0 Multiplier on artificial density 
correction term, to increase damping. 
Maximum, 2.0. Intended to reduce 
overshoot on midblade shocks. 

0. 2 To improve stability for multiple-grid 

cases, SUPREL is reduced by an amount 
equal to CII for each successive grid 
until SUPREL = 0.5. 

10 Number of iterations between intermediate 

printouts of residuals and Mach number 
and between updates of the densities 
used in equation (21) . 

FALSE =TRUE. to list flow quantities at all 
grid points on the last mesh. 

FALSE =.TRUE. to activate stream-channel 

thickness and/or radius variations. 

1 Logical position of desired stream- 

surface data on the stream-channel 
file. See the section Bulk data for 
flow solution runs. 

1. Spanwise radius at upstream boundary 
divided by aerodynamic chord. Radius 
effects are activated if RINF ^ 1. With 
the current version of the code, the 
following cases are allowed: 


QUASI 3 

RINF 

Results 

.F. 

1. 

Planar two-dimensional 

.T. 

^1. 

Thickness on file 
Radius on file 

.T. 

1. 

Thickness on file 
Constant radius 


0. Thickness of a wake, as a fraction of 

aerodynamic chord. Used only in 
calculation of downstream Mach number 
to simulate a blockage. 
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MINF2 


Real 


OMEGA Real 

VAXIAL Real 


FLOCO Real 


10.0 


0 . 


999. 


999. 


Downstream Mach number. If any value is 
input other than 10., it is used as 
Mout* Otherwise, Mont is calculated by 
using equation (27). MINF2 is available 
as another means of introducing blockage 
(see WAKE) or to compare a run with 
other analyses. Also, because of the 
behavior of the function F in equation 
(27) near sonic speeds, it is often 
necessary to input the downstream Mach 
number here and to avoid use of equation 
(27) if Mont "is known to be near 1.0. 

For rotors, wheel rotation is in radians 
per second. If blade row is moving down 
along y axis, OMEGA < 0. 

Must be input if OMEGA ^ 0. VAXIAL is 
the upstream velocity component in the 
streamwise (x) direction divided by the 
aerodynamic chord. VAXIAL = V /C. 

/ V 1 

Units are ( seconds )“■*■ . 


For rotors, wheel speed can be input by 
using OMEGA and VAXIAL or as FLOCO, the 
flow coefficient at the upstream 
boundary. FLOCO = /U, where U is 

wheel speed. If OMEcA'^is used for input, 
FLOCO should be left as 999. It will be 
calculated as 


FLOCO = 


VAXIAL 

OMEGA • RINF 


Bulk Data Input 

Bulk data for geometry runs. - The format for all variables listed here 
is namelist form. The namelist is INSTUF. 


Name Type Default Description 

H2 Complex NONE Table of all blade definition points. 

The real part is the x coordinate and 
the imaginary part is the y coordinate. 
The table should begin at the point of 
maximum x value (near the trailing 
edge) and proceed clockwise back around 
to the first point, which is repeated. 

The blade must be at stagger angle and 
the origin near the point of minimum 
X. For electrostatic analog grids the 
stagger angle must be positive (i.e.. 
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the leading edge is low, and the 
trailing edge is high). The maximum 
number of points in H2 is 99 for the 
interpolation scheme or 63 for the 
electrostatic analog. 

BUG2 Logical FALSE = .TRUE, for a more detailed output list 

of the geometry generation run, 
including second derivatives at grid 
points on the body. 

B ulk data for flow solution runs. - The mesh point storage file (unit 
23) is a required input only if a mesh is not generated in the same run as 
the flow solution. The stream-channel data file (unit 2) is required only 
if quasi-three-dimensional effects are to be calculated. 

Mesh point storage file: If the mesh point file is needed, the values 

of NOZES to be used must be negative. The file is written by a geometry 
generation run as unit 13. The first record is 215, with the first position 
being the value of NOZES for which the grid was generated and the second 
position being the value of MS used. All other records are formatted 6E13.5 
and contain the x values at all points in the mesh, followed by all y 
values. The sequence for both is to begin on the periodic boundary at the 
line going to upstream infinity (nominal) and then to proceed clockwise 
around each ellipse-like curve, repeating the first point, moving in to the 
next ellipse, and continuing until the blade surface has been listed in the 
same way. The orientation of the x,y points is with the aerodynamic chord 
along the x axis and the origin at midchord. The statements that create 
this file in subroutine WRAPUP have the form: 

WRITE (13,1) NOZ, M, KN, NED 

WRITE (13,2) ((X(I,J), I = 1, NZZ), J = 1,M), ((Y(I,J),I = 1, NZZ), 
d = 1,M), (ETA(I), I = 1, NZZ), (ZETA(J), J = 1,M) 

1 FORMAT (415) 

2 FORMAT (6E13.5) 

Stream-channel data file: If the stream-channel data file is required, 

QUASI3 must be .TRUE. The file must contain a table of stream-tube 
thicknesses, radial positions, and corresponding x values along the stream 
surface. The format is identical to that shown on page 18 of reference 2. 
However, not all of these data are necessary for QSQNIC. The following list 
and figure 8 show the minimum requirements for the data needed by QSONIC. 

If more than one such set of data exists on the file, NSTRM is used to 
select the desired set. (See the section TEST CASE.) 


Record 

Column 

Name 

Type 

Default 

Description 

number 

number 





1 

BLANK 





2 

BLANK 





3 

21-30 

CHO 

REAL 

NONE 

Aerodynamic chord 
multiplied by the cosine 
of LAMDAO 

4 

BLANK 
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5 36-40 NRSP INTEGER NONE Total number of data 

points in each of the 
tables of thickness, 
radial position, and x 
locations. If NRSP = 2, 
a linear distribution is 
obtained between the end 
points given. 

6 BLANK 

7 

8 

9 

10 

11 


12 

1-80 RM 

(8F10.5) 

REAL 

NONE 

Array of corresponding 
locations for thickness 

to 



and radius data values. 

conclusion 




For a blade sitting at 
the proper stagger angle 
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X = 0 at the blade leading 
edge. Units can be any 
consistent length scale 
common to CHO, RM, RMSP, 
and BESP. 



1-80 

(8F10.5) 

RMSP 

REAL 

NONE 

Spanwise radial positions 
of stream surface at given 
X locations. 

LAST 

1-80 

(8F10.5) 

BLANK 

BESP 

REAL 

NONE 

Array of stream-tube 
thickness values at the 
given x locations. 


Restart input file: This is required if RESTAR = .TRUE. This file was 

written by the code during the interrupted run. It contains the field of 
potential values in 6E13.5 format and is updated whenever a new grid is 
begun (if more than one was requested with MESHl ^ MESHN) and also at 
intervals of the iteration count divisible by the input parameter IT. 

Caution: Restart capabilities have not been fully tested in version 1.5 

of the code. 


Output 


All output produced by this program falls into the following categories: 

(1) Geometry generation run output 

(a) Printed listing (unit 6) 

(b) Mesh points file (unit 13) 

(2) Flow solution run output 

(a) Printed listing (unit 6) 

(b) Restart save file (unit 8) 

(c) Plot data save file (unit 13) 

Refer to the test case in the next section for examples of the printed list 
outputs. 

Geometry Generation Run Outputs 

Printed listing. - The output list on unit 6 will contain the following 
items: 


(1) Title of case, as entered on input 

(2) The phrase CASE CONTROL INPUT ECHO, followed by a namelist output 
of the input data as entered with NAMELIST PARAMS. All variable names are 
identical to those described in the INPUT section. 

(3) For the interpolation scheme the phrase INTERPOLATION SCHEME GRID 
GENERATED, followed by any conditional messages as described in the section 
on error messages 

(4) For the electrostatic analog the phrase ELECTROSTATIC ANALOG GRID 
GENERATED, followed by 

(a) The phrase NEGRID = , followed by the value of MS used for this 
grid 

(b) The phrase NZGRID, ODD = , followed by the value of NOZES used 
for this grid. 
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(c) The phrase IMAXY = , followed by the index of the point with 
the maximum y value 

(d) The phrase SURFACE COORDINATES, followed by a two-column list 
(x, then y) of the blade coordinates as obtained from NAMELIST 
INSTUF 

(e) The phrase NUMBER OF BODY POINTS =,, followed by the total 
number input from INSTUF 

(f) The phrase PITCH = , followed by input value of S divided by 
CHORD 

(g) The phrase CHO = , followed by the value of CHORD * COS(LAMDAO) 

(h) The phrase STAGGER = , followed by the value of LAMDAO in 
degrees 

(i) The phrase ITER COUNT FOR INTERNAL GRID POINTS = , followed by 
the number of iterations required to obtain all points between 
the periodic boundary and the blade 

(j) The phrase ERR = , followed by the average residual for this 
iteration 

(5) If BUG2 = .TRUE., these additional outputs and phrases will be 
listed for either grid generator: 

(a) NOZ = , followed by the value of NOZES used 

(b) M = , followed by the value of MS used 

(c) KN = , followed by number of blade definition points on upper 
surface of blade after going through the smoothing routine 
Smoothing done only for Interpolation scheme grids) 

(d) NED = , followed by the total blade definition points after 
smoothing 

(e) A list of smoothed blade definition points along with 
calculated second derivatives. For electrostatic analog grids, 
blade orientation is the same as described under parameter H2 
in the input section. For interpolation scheme grids, it is 
the same as described in item (5h). 

(f) All X grid points are listed as they appear on the mesh points 
storage file 

( 9 ) All y grid points are then listed in a similar manner 

(h) Final calculated grid points that lie on the blade are listed 
along with calculated second derivatives. The orientation of 
the blade (interpolation scheme) is with the chord horizontal 
and the origin at midchord. For electrostatic analog grids, 
orientation is the same as on input. 

(6) Final calculated grid points that lie on the blade are listed. 

Orientation is with the chord at stagger and the origin at midchord. 

Mesh points file. - The mesh points file is written on unit 13 and 

contains one record (215) with the values of NOZES and MS used in generating 
the grid, followed by all x values, then all y values, in 6E13.5 format. No 
user intervention is required to use the file as bulk input for a future 

flow solution run (unit 23). It is also useful as plot input for 

preliminary grid views. 

Flow Solution Run Output 

Printed listing. - The printed listing is similar to grid generation 
runs for the first two items: 
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(1) Title of case 

(2) The phrase CASE CONTROL INPUT ECHO, followed by namelist list 

(3) For cases with rotation (OMEGA > 0. or FLOCO ^ 999.) the phrase 
ROTATION EFFECTS INCLUDED, ABSOLUTE AND RELATIVE TRIANGLES HAVE -, 
followed by a NAMELIST output of the following velocity triangle 
parameters: 

MINF Input value of relative Mach number (relative to the frame 
fixed to blade) 

MINF2 Relative Mach number downstream as calculated or input 

QINF2 Relative speed downstream, q relo^t/Q >^®^in ^ ’'®^in ~ 

where q = (u^ + v^)^/^ 

ABETAl Absolute inlet flow angle measured from throughflow (x) 
direction to actual flow direction, where clockwise is 
negative 

ABETA2 Absolute outlet flow angle at downstream boundary 

AQINFl Upstream speed ratio, absolute to relative, q abs-j^/Q '"e^in 

AQINF2 Downstream speed ratio, q abSout/9 '^®^in 

AMINF Absolute upstream Mach number (subtracting wheel speed from 
MINF) 

AEM Absolute downstream Mach number 

ROTATN Wheel speed as calculated from FLOCO and BETAl normalized by 
assuming q relin = 1. ROTATN = cos (BETAl + LAMDAO) /FLOCO. 

For all flow cases: 

(4) The phrase CALCULATED FLOW PARAMETERS, followed by a Namelist 
output of these calculated downstream conditions: 

CIRCLN Exact circulation for given inlet and outlet flow angles 

MINF2 Relative Mach number downstream as calculated or input 

B20B1 Ratio of downstream channel thickness to upstream thickness, 
as input or as obtained from stream-channel data file 

R20R1 Ratio of downstream spanwise radial position of stream- 

channel to upstream position, as input or as obtained from 
stream-channel data file 

RINF Spanwise radius at upstream boundary divided by aerodynamic 

chord, as input 
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RATLE 


Calculated radius at leading edge of blade based on radius 
distribution from stream-channel data file or the linear 
distribution 


For each grid selected for relaxation, 


(5) A table of intermediate results is printed at intervals whenever 
the iteration count is divisible by IT. The table consists of the 
following: 

N 

AVERAGE PHI CORRECTN ^ 

1 



RESIDUAL AT T*E* 
MAX MACH ON BLADE 


Net mass flow out of finite volume element 
at trailing edge 

Mmax blade surface', not including points 
in the first and last 6 percent of blade chord 


CALCULATED CIRCULATION C^ = (<P]^st “ ‘•’l)on blade 

line 


RELATIVE CIRC ERROR 
SUBSONIC RELAX FAC. 
SUPERSONIC RELAX FAC. 


(Cc - /^exact> where - CIRCLN 

OVEREL, as input, subject to delay by NOWREL 

SUPREL, as input, subject to reduction by Cl I 
and delay by NOTYET 


(6) For each grid selected for relaxation, a table of final flow 
calculation on blade surface, including X/CX, Y/CX, X, Y, MACH, STATIC 
PRES COEF, PHI, XVEL, YVEL, and DENSITY, where CX = CHORD * (COS 
LAMDAO), PHI denotes the potential field values, XVEL and YVEL are 
velocity components normalized to u^p + v^p = 1, and density 

is normalized to p-jp = 1* 

(7) For each grid selected for relaxation, after all grids have been 
solved, and if ALLOUT = .TRUE., a table of flow variables at all grid 
points is printed, beginning with the phrase CALCULATED FLOW FIELD AT 
ALL GRID POINTS. All variables are the same as in the previous table 
(item 6) with the addition of FLOW ANGLE = tan“^ YVEL/XVEL in degrees. 


Restart save file. - See section Bulk data for fl ow s ol utio n s runs for a 
description of this file, written on unit 8. 


Plot data save file. - At the completion of the table of the flow field 
at all grid points, or immediately after the last grid is complete if ALLOUT 
= .FALSE., the following information is recorded on unit 13: 
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Record 

number 


Format 


Data 


1 215 N0Z,M values of NOZES and MS on final grid used. 

Total number of grid points = (2 * NOZ - 1) * (M) = 
NOO * M = NTOT 

NTOT gnio c ({X(I,0),I = 1,N00),J = 1,M), throughflow coordinates. 
"6 * Blade is at stagger. Origin is at midchord. 

((Y(I,J),I = 1,N00),J = 1,M), tangential coordinates 


((XVEL(I,d),I = N00),0 = 1,M), throughflow velocity 
components 

((YVEL(I,J),I = NOO),J = 1,M), tangential velocity 
components 

((MACH(I,J),I = N00),J = 1,M), Mach numbers 


( (PRES(I, J) ,I = N00),J = 1,M), static pressure 
coefficients 

If a flow solution and geometry generation are combined in one run, this 
plot data will overwrite the mesh points (also unit 13) written by the grid 
generator. The graphic output shown with the sample test case was generated 
by a local graphics package at NASA Lewis that accessed a PLOT data file 
such as this. 

Error and Condition Messages 

Module Message phrase Description 

XYCALC SGEN FAILED. IBAD= In attempting to smooth a table of input 

data, the module that calculates the value 
of the distance along the curves detected a 
large change in slope between adjacent point 
pairs at the table index number of IBAD. 
Remedy; If the output grid points are 
unsatisfactory, check the data going into 
XYCALC for any points out of sequence, 
repeated points, all-zero data, etc. 

XYCALC SPGEN UNABLE TO The module that calculates new points on 

COMPLETE SEGMENT the curve to smooth the data could not 

satisfy the point spacing criteria listed in 
the SPGEN subroutine. Remedy: Improve the 

input data (i.e., add more points if 
necessary) so as to avoid any large changes 
in slope (>90“) between adjacent points. 
Also, check the data for discrepancies as 
noted for the previous message. 


o NTOT 

* 2 -5— 

, NTOT 

* 3 

.4^ 
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The following four messages refer to line .numbers While a mesh is 

being generated, radiating lines are numbered beginning with the line to 
upstream infinity and proceeding clockwise around the blade. 


SURFUP ITERATION PROCEEDS 
WITH SMOOTHER 
SLOPE DATA FOR 
LINE 


After 20 attempts to find the intersection 
of the blade surface and radiating line 
the table of x and y values versus n 
(surrounding contours) on line 'I' is 
smoothed by using subroutine XYCALC, and the 
iteration is attempted again. 


SURFUP ITERATION FAILED 
FOR LINE 'I'. 
PROCEEDING TO EDGE 
REGION TECHNIQUE 


When the total iteration attempts mentioned 
above reach 50 or they become divergent, the 
solution of this line is put aside until all 
others are completed. At that time the 
intersection point is found by interpolating 
on the table of intersection points formed 
by all the successful iterations. This is 
referred to in the code as the "edge region 
technique," because this is normally 
required in the leading- and trailing-edge 
regions of sharp blades. 


SURFUP ITERATION NOT 
ATTEMPTED FOR 
LINE ■!'. BODY X 
TOO CLOSE TO LE 
OR TE: 'XB XP' 

PROCEED TO EDGE 
TECH." 


(See explanation of CHOP parameter in the 
section on geometry parameters.) If X = XB 
(where XB is the intersection of radiating 
line with the X axis) is larger in 
absolute value than the quantity CH0P*XP 
(where XP is equal to the half-chord length 
of the flat plate), this line is assumed to 
be near a sharp tip region on the blade and 
the normal iteration procedure is bypassed. 
This message should be expected if CHOP < 
0.95. 


SURFUP GRID LINE NO. 

HAS DEGENERATED 
TO POINT. TRY 
THETL,THETT = 
LE,TE CAMBER 
ANGLES. 


For certain high-camber blades with small tip 
radii, the stretching of the grid to include 
the blade can degenerate if the camber 
information is not available at the leading 
and trailing edges. See THETL and THETT 
parameter definitions in the section on 
geometry inputs. 


All other messages are self-explanatory or not critical to code success. 


TEST CASE 

The execution times for the following sample case, run on an IBM 3033, were 

(1) Interpolation scheme grid generation, 40 seconds for 2500 grid 
points or 

(2) Electrostatic grid generation, 30 seconds and 

(3) Quasi-three-dimensional flow solution on three consecutively finer 
meshes, 260 seconds 
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The blade geometry 1s the hub section of the Sanz supercritical compressor 
stator designed for NASA Lewis by using the methods of reference 6. A good 
compari son between QSONIC and the original hodograph design solution is 
presented in reference 3 for a planar two-dimensional flow. 


Input - Case Control Parameters for Interpolat ion Scheme Grid Generation 


SANZ SUPERCRITICAL STATOR HUB SECTION, GRID GENERATION RUN 
CPARAMS S= . 86 3 , NSTRM=3 , RINF= 3 . 72<t , 

MESH1 = ^,MESHN = 4,QUASI3=:. FALSE. , 

MS=3,7,13,25,NOZES = -9,-13,-25,<i9, 

NOFLOM=.TRUE. , 

TOLS=l .E04,-2.E-05, 

7, E-05, 6. E-05,-. 000001,2*0. , 

BETA1=25. 1 ,NOMREL=20 ,NOTYET=10, 

BETA2=-15.6,LAMDA0=19.5, 

OVEREL= 1.6, CII= . , MINF= .801, 

IT=3 0 , ALLOUT= . TRUE . , SUPREL=2 . 0 , SMOOTH= . TRUE . 

RLE=. 007,RTE=.029,KN=62,NED=88,CHORD=1 .01^6 


Bulk Data Input for Geometry Generation 


8INSTUF 

Hl= 100»<(0 . 0, 0.0), 

H2= <0.9^980, 0.3523999), (0.9^82999, 0.3^2^0), (0.93^2999, 0.325^0), 

( .9182999,0.321^0), 

(0.9072999, 0.321^999), (0.8672999, 0.32190), (0.803^999, 0.32310), 

( .6818999,0.31^4), 

(0.5879999, 0.29590), (0.5217999, 0.27630), (0.4682999, 0.2560), 

( .4284999, .23840), 

(0.3942999, 0.22160), (0.3675999, 0.2072999), (0.3058999, 0.17120), 
(.2519999, .13680), 

(0.2052999, 0.10620), (0.1571999, 0 . 7459998E-01 ) , (0.1080999, 0.4380E-01), 

(0 .8219993E-01 , 0 .2860E-01 ) , (0 .6140E-01, ' 0.1680E-01 ) , 
(0.03989999,0.006399997), 

(0.2889999E-01, 0 . 1 1 9999 9E-02 ) , ( 0 , 1789999E-0 1 , -0.250E-02), 

( . 01259999,-. 0036) , 

(0.7299997E-02, -0.360E-02), (0.0, 0.0), (-0 . 1000017E-03, 0 . 71 99995E-02 ) , 

( 0. 3999 9 97 E- 02, 0.2030E-01), (0.160E-01, 0.4080E-01), 

( .03529999, .06469995), 

(0.6869996E-01, 0.1005999), (0.1130999, 0.14380), (0.1557999, 0.18190), 

(0.1867599, 0.207410), (0.20840, 0.22420), (0.225390, 0.236780), 

( .2392299, .246570), 

(0.2507299, 0.2543899), (0.2607099, 0.260950), (0.26950, 0.2665499), 
(.2773799, .271430), 

(0.2844899, 0.2757199), (0.2909699, 0.279530), (0.2969099, 0.282940), 

( .3023999, .286020), 

(0.3074799, 0.288810), (0.3122399, 0.291370), (0.316680, 0.293710), 
(.3208599, .2958699), 

(0.3247699, 0.297850), (0.3285099, 0.299720), (0.3320799, 0.301460), 
(.335510, .3031099), 

(0.33880, 0.304660), (0.3419899, 0.3061399), (0.3450699, 0.307550), 

( . 3480799, . 30890) , 

(0.3510199, 0.31020), (0.3539199, 0.311460), (0.356750, 0.312660), 

(. 3595399, .313840), 

(0.3622699, 0.3149599), (0.364960, 0.316050), (0.3675599, 0.317090), 
(.3701199, .318090), 

(0.37260, 0.3190399), (0.375010, 0.319950), (0.3773299, 0.320810), 
(.3795499, .321620), 

(0.3816999, 0.322390), (0.3838699, 0.323140), (0.3861299, 0.323920), 
(.4001899, .328540), 

(0.40740, 0.33060), (0.4604999, 0.34350), (0.5118999, 0.35270), 

( .5662999, .36020) , 

(0.6176999, 0.36580), (0.6926999, 0.37220), (0.7486999, 0.37590), 

( .8201999, .3795), 

(0.8772999, 0.38140), (0.9072999, 0.38230), (0.9172999, 0.38250), 

( ,93330. .3794), 

(0.94730, 0.36440), (0.94980, 0.3523999), 12»(0.0, 0.0), 

BUG2= F, 

tEND 
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Printed Output List for Geometry Generation 


QSOHIC VERSION 1.5 SANZ SUPERCRITICAL STATOR HUB SECTION, GEOMETRY GEN 

CASE CONTROL INPUT ECHO 

ePARAMS 
NOFLOW= T 
RE5TAR= F 
PEMESH= 4 
MS= 3, 7, 13, 25, 

HOZES= -9, -13, -25, ^9, 6*0 

MESH1= 4 

ME5!!N= 

LAMDA0= 19.50 
CHORD= 1.01<V60 
S= 0.B630 
NED^ 8S 
KK= 62 

RLE= 0.6999999E-02 
RTE= 0.290E-01 
THETL= 0.0 
THETT= 0.0 
CAMPER= 6 
STABAC= 0.9990 
CHOP= 0.9899999 
SMOOTH= T 
LEONLY= F 
SLP1= 1.0 
SLP2= 2.0 
SLP3= 1.0 
SLP^= 1.0 
MINF= 0.8010 
BETA1= 25.09999 
BETA2= -15.60 
GAnriA= 1.^0 

TOLS= 10000.0, -0.1999999E-0<i, 0 . 6 999999E-0 A , 0 . 5999999E-0 A , -0 . 999999<»E-06 
5*0.0 

OVEREL= 1.599999 
UNDERL= 1.0 
SUPREL= 2.0 
NOWREL= 20 
NOTYET= 10 
TE€ARD= 2.0 
DAMP= 1.0 
CII= 0.^0 
1T= 30 
ALLOUT= T 
QUASI3= F 
NSTPri= 3 
RIHF= Z.72^0 
UAKE= 0.0 
OMEGA= 0.0 
VAXIAL= 0.0 
FLOCO= 999.0 
riINr2= 10.0 
BIGLAK= F 


INTERPOLATION SCHEKE GRID GENERATED 


SCEN FAILED. IBAD> 10 

ITERATION PROCEEDS WITH SMOOTHER SLOPE DATA FOR LINE 
ITERATION PROCEEDS UITH SMOOTHER SLOPE DATA FOR LINE 4<i 

ITERATION FAILD FOR LINE 4A PROCEEDING TO EDGE REGION TECHNIQUE 


ITERATION 

NOT 

ATTEMPTED 

FOR 

LINE 

45, 

BODY 

X 

TOO 

CLOSE 

TO 

LE 

OR 

TE: 

0.488165E 

00 

0.491130E 

00 

PROCEED 

TO 

EDGE 

TECH. 

ITERATION 

NOT 

ATTEMPTED 

FOR 

LINE 

46, 

BODY 

X 

TOO 

CLOSE 

TO 

LE 

OR 

TE: 

0.489080E 

00 

0.4ni30£ 

00 

PROCEED 

TO 

EDGE 

TECH. 

ITERATION 

HOT 

ATTEMPTED 

FOR 

LINE 

47, 

BODY 

X 

TOO 

CLOSE 

TO 

LE 

OR 

TE: 

8.487084E 

00 

0.49I130C 

• 0 

PROCEED 

TO 

EDGE 

TECH. 

ITERATION 

NOT 

ATTEMPTED 

FOR 

LINE 

93, 

BODY 

X 

TOO 

CLOSE 

TO 

LE 

OR 

TE: 

-0.488160E 

00 

0.491130E 

00 

PROCEED 

TO 

EDGE 

TECH. 

ITERATION 

HOT 

ATTEMPTED 

FOR 

LINE 

94, 

BODY 

X 

TOO 

CLOSE 

TO 

LE 

OR 

TE: 

•0.489081E 

00 

0.491130E 

00 

PROCEED 

TO 

EDGE 

TECH. 

ITERATION 

HOT 

ATTEMPTED 

FOR 

LINE 

95, 

BODY 

X 

TOO 

CLOSE 

TO 

LE 

OR 

TE: 

-0.487099E 

00 

O.49m0E 

00 

PROCEED 

TO 

EDGE 

TECH. 


SGEH FAILED. IBAD= 21 
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GRID OUTPUT POINTS ON BLADE SURPACE, CHORD AT STAGGER ANGLE 
I X y 


0 

0.95668E 00 

0.15691E 00 

1 

0.99927E 00 

0.15395E 00 

2 

0.93797E 00 

0.15170E 00 

3 

0.92350E 00 

0.15180E 00 


0.90665E 00 

0.15197E 00 

5 

0.38739E 00 

0.15221E 00 

6 

0.36596E 00 

0.15272E 00 

7 

0.39269E 00 

0.15328E 00 

8 

0.31783E 00 

0.15326E 00 

9 

0.29187E 00 

0.15299E 00 

10 

0.26999E 00 

0.15099E 00 

11 

0.23798E 00 

0.19860E 00 

12 

0.20961E 00 

0.19530E 00 

13 

0.18160E 00 

0.19100E 00 

U 

0.15360E 00 

0.13576E 00 

15 

0.12576E 00 

0.12961E 00 

16 

0.98172E-01 

0.12256E 00 

17 

0.70919E-01 

0.11969E 00 

18 

0.99072E-01 

0.10585E 00 

19 

0.17699E-01 

0.96168E-01 

20 

-0.82057E-02 

0.85733E-01 

21 

-0.33653E-01 

0.79616E-01 

22 

-0.58602E-01 

0.62837E-01 

23 

-0.83080E-01 

0.50958E-01 

29 

-0.10701E 00 

0.37993E-01 

25 

-0.13050E 00 

0.29196E-01 

26 

-0.15353E 00 

0.10958E-01 

27 

-0. 17608E 00 

-0.39693E-02 

28 

-0.19820E 00 

-0.17969E-01 

29 

-0.21983E 00 

-0.31533E-01 


30 -0.2<tl09E 00 -0.45«52E-01 

31 -0.26181E 00 -0.59060E-01 


32 

-0,28200E 

00 

-0.72396E 

-01 

33 

-0.30169E 

00 

-0.85332E 

-01 

39 

-0.32071E 

00 

-0.97631E 

-01 

35 

-0.33911E 

00 

-0.10929E 

00 

36 

-0.35678E 

00 

-0.12017E 

00 

37 

-0.37360E 

00 

-0.13018E 

00 

38 

-0.38993E 

00 

-0.13993E 

00 

39 

-0.90911E 

00 

-0.19766E 

00 

90 

-0.91773E 

00 

-0.15998E 

00 

91 

-0.92976E 

00 

-0.16019E 

00 

92 

“0.99028E 

00 

-0.16992E 

00 

93 

-0.99929E 

00 

-0.16776E 

00 

99 

-0.95658E 

00 

•0.16885E 

00 

95 

-0.96225E 

00 

-0.16779E 

00 

96 

-0.96695E 

00 

-0.16957E 

00 

97 

-0.96638E 

00 

-0. 15896E 

00 

98 

-0.96950E 

00 

“0.15070E 

00 

99 

-0.96095E 

00 

-0.19111E 

00 

50 

-0.95376E 

00 

-0.12980E 

00 

51 

-0.99955E 

00 

-0,11689E 

00 

52 

-0.93276E 

00 

-0.10285E 

00 

53 

“0,91890E 

00 

-0.87286E- 

'01 

56 

-0 .90 J09E 

00 

•O.rOftOSE-fll 

55 

-0 . 38557E 

00 

-0.52916E-01 

56 

-0.36655E 

00 

•0.39965E- 

'01 

57 

-0.39618E 

00 

-0.15981E- 

-01 

58 

“0.32971E 

00 

0.37890E- 

-02 

59 

-0.30223E 

00 

0.22980E- 

-01 

60 

-0.27893E 

00 

0.91877E- 

-01 

61 

-0.25996E 

00 

0.60262E- 

•01 

62 

-0.23093E 

00 

0.77926E- 

■01 

63 

-0.20591E 

00 

0.99578E- 

■01 

69 

-0.17998E 

00 

0.11009E 

00 

65 

-0.15921E 

00 

0.12932E 

00 

66 

-0.12812E 

00 

0.13713E 

00 

67 

-0.10166E 

00 

0.19815E 

00 

68 

-0,79903E- 

■01 

0.15798E 

00 

69 

-0.97999E- 

■01 

0.16538E 

00 

70 

-0.20778E- 

■01 

0.17171E 

00 

71 

0.63922E- 

■02 

0.17716E 

00 

72 

0.33521E- 

■01 

0.18183E 

00 

73 

0.60565E- 

■01 

0.18581E 

00 

79 

0.87928E- 

■01 

0.18939E 

00 

75 

0.11911E 

00 

0.19297E 

00 

76 

0.19052E 

00 

0.19529E 

00 

77 

0.16669E 

00 

0.19770E 

00 

78 

0.19239E 

00 

0.19989E 

00 

79 

0.21772E 

00 

0.20189E 

00 

80 

0.29257E 

00 

0.20358E 

00 

81 

0.26689E 

00 

0.20510E 

00 

82 

0.29097E 

00 

0.20697E 

00 

83 

0.31337E 

00 

0.20765E 

00 

89 

0.33599E 

00 

0.20867E 

00 

85 

0.35660E 

00 

0.20999E 

00 

86 

0.37672E 

00 

0.21015E 

00 

87 

0,39582E 

00 

0.21079E 

00 

88 

0.91396E 

00 

0.21131E 

00 

89 

0.93138E 

00 

0.21187E 

00 


I 
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90 

0.9^81^E 

00 

0 .21073E 

00 

91 

0.46156E 

00 

0.20281E 

00 

92 

0.46797E 

00 

0.19276E 

00 

93 

0.^6992E 

00 

0.18158E 

00 

9<* 

0,<f6799E 

00 

0.17112E 

00 

95 

0.<»6272E 

00 

0.162^8E 

00 

96 

0.^5668E 

00 

0.15691E 

00 

97 

0.4^927E 

00 

0,153^5E 

00 

98 

0.43797E 

00 

0.15170E 

00 


Case Control Parameters for Flow Solution 


SANZ SUPERCRITICAL STATOR HUB SECTIOH, FLOW SOLUTIOH 
CPARAnS S» . 86 3 , NSTRM=3 , RINF=3 .726* 
t 1 ESHl = 2,nESHN=6,QUASI3=.TRUE. , 
«S=3,7,13.25,NOZES=-9.*13.-25,-69, 

TOLS=1.E09,-2.E-05, 

7,E-05.6.E-05,-.00 0001.2»0. , 

BETA1=25.1,NOWREL=20,NOTTET=10. 

BETA2=-15.6,LAMDA0=19.5, 

OVEREL= 1 . 6 , CII= . <t , MINF= . 8 0 1 , 

IT= 3 0 , ALLOUT= . TRUE . . SUPREL=2 . 0 , SMOOTH* . TRUE . 

RLE*. 0 0 7, RTE*. 029, KN*62.NED=88, CHORD®!. 0196 
£END 


Bulk Data Input for Flow Solutions - Stream-Channel Data File 


"Streamline number 3" is selected here because NSTRM equals 3 from the 
case control input. 


STREAMLINE NUMBER 1 -- STREAM FUNCTION = 0.0000 


l.<t00 00 

1716.4800 

581.9998 

0.00278050 

1.00061536 


0.000 

0.00000 

290.7671 

-4.9297 

0.119776 

0.1001226 





1.0000 
21 61 

0.0010 

81 20 

28 28 

1 1 

1 



0.0013696 

0.0038565 

54.3510 

1.8563 

13.0000 




0.000256 

0.007862 

0.018306 

0.028770 

0.039246 

0.049726 

0.060202 

0.070664 

0.0S1105 

0.091519 

0.101903 

0.112255 

0.115792 




0.0017591 

0.0239841 

0.0491248 

0.0677468 

0.0829612 

0.0930837 

0.0986944 

0.1022954 

0.10^7776 

0.1063557 

0.1075357 

0.1080631 

0.1082771 




0.0013696 

0.0038565 

16.9621 

-2.6960 

13.0000 




0 .00 1 768 

0.007862 

0.018306 

0.028770 

0.039246 

0.049726 

0.060202 

0.070664 

0.081105 

0.091519 

0.101903 

0.112255 

0.115736 




■0.0028870 

0.0023670 

0.0170644 

0.0341531 

0.0495386 

0.0626676 

0.0736161 

0.0817172 

0.0875790 

0.0911617 

0.0924746 

0.0922745 

0.0919729 




-0.150957 

-0.118756 

-0.086551 

-0.054317 

-0.043995 

-0.033661 

-0.023312 

-0.012945 

-0.002555 

0.007863 

0.018306 

0.028770 

0.039246 

0.049726 

0.060202 

0.070664 

0.081105 

0.091519 

0.101903 

0.112255 

0.122580 

0.132888 

0.143187 

0.163188 

0.183188 

0.203190 

0.223192 

0.243194 





0.<i«»5770 

0.445983 

0.446553 

0.448051 

0.448800 

0.449698 

0.450754 

0.451980 

0.<»5338<» 

0.454976 

0.456732 

0.458606 

0.460548 

0.462512 

0.464450 

0.466314 

0.<»68056 

0.469629 

0.470984 

0.472074 

0.472866 

0.473389 

0.473691 

0.473822 

0.4736^1 

0.473389 

0.473108 

0.472820 





0.0060057 

0.0068080 

0.0073119 

0.0073806 

0.0073717 

0.0073965 

0.0073621 

0.0074961 

0.0069645 

0.0071813 

0.0064250 

0.0060114 

0.0060900 

0.0058698 

0.0055469 

0.0053046 

0.0051044 

0.0051718 

0.0050849 

0.0049574 

0.0046936 

0.0049132 

0.0048774 

0.0050551 

0.0051724 

0.0051911 

0.0052327 

0.0052886 





0.74436 

0.71080 

0.69488 

0.69051 

0.68945 

0.68729 

0.68636 

0.68103 

0.69274 

0.74329 

0.71332 

0.68195 

0.56593 

0.49017 

0.46927 

0.46599 

0.47064 

0.45210 

0.45687 

0.47149 

0.46048 

0.43514 

0.43864 

0.42024 

0.40927 

0.40781 

0.40432 

0.39965 





0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00143 

0.00333 

0.00524 

0.00714 

0.00905 

0.01096 

0.01286 

0.01476 

0.01666 

0.01855 

0.02043 

0.02180 

0.02180 

0.02180 

0.02180 

0.D2180 

0.02180 

0.02180 

0.02180 





I 0 

0 2 

2 2 

3 





STREAMLINE NUMBER 2 -- 

STREAM FUNCTION * 0 

.0290 



1.40000 

1716.4800 

581.7185 

0.00300310 

.00061536 


0.000 

0.00000 

299.9548 

15.3325 

0.120048 

0.0929822 





1.0000 
21 61 

0.0010 

81 20 

28 28 

1 1 

1 



0.0015167 

0.0038388 

53.7777 

2.3267 

13.0000 




0.000293 

0.005101 

0.015332 

0.025679 

0.036085 

0.046562 

0.057124 

0.067719 

0.078345 

0.088996 

0.099678 

0.110391 

0.116051 




0.0019039 

0.0153618 

0.0396086 

0.0585241 

0.0737625 

0.0847238 

0.0913482 

0.0952344 
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0.0977184 0.0990426 0.0999755 0.1004877 0.1008875 

0.0015167 0.0038388 14.8797 *1.8425 13.0000 

0.001904 0.005101 0.015332 0.025679 0.036083 0.046562 0.057124 0.067719 

0.078345 0.088996 0.099678 0.110391 0.116083 

*0.0031139*0.0008983 0.0115014 0.0267384 0.0413744 0.0544321 0.0656631 0.0744360 
0.0807123 0.0841923 0.0855135 0.0854021 0.0850744 

-0.151271 -0.119060 -0.087167 -0.055544 -0.045453 -0.035364 -0.025275 -0.015178 

-0.005052 0.005101 0.015332 0.025679 0.036083 0.046562 0.057124 0.067719 

0.078345 0.088996 0.099678 0.110391 0.121096 0.131718 0.142263 0.162577 

0.182676 0.202704 0.222718 0.242718 

0.460865 0.462438 0.463745 0.465257 0.465955 0.466842 0.467804 0.469145 

0.469790 0.471835 0.472950 0.473955 0.475829 0.477420 0.478547 0.479806 

0.481111 0.482881 0.484189 0.485120 0.485247 0.486319 0.486553 0.487007 

0.487043 0.486819 0.486621 0.486447 

0.0046231 0.0048699 0.0050106 0.0050022 0.0049897 0.0049849 0.0049670 0.0049825 

0.0048608 0.0050024 0.0049863 0.0048432 0.0049296 0.0047053 0.0044550 0.0042884 

0.0041698 0.0042085 0.0042049 0.0041758 0.0039849 0.0041184 0.0041017 0.0041611 


0.0042039 

0.0042081 

0.0042270 

0.0042528 





0.80177 

0.77472 

0.76107 

0.75856 

0.75805 

0.75658 

0.75593 

0.75210 

0.75994 

0.79923 

0.79827 

0.78752 

0.69086 

0.61466 

0.57647 

0.55905 

0.55276 

0.53098 

0.52659 

0.53366 

0.52112 

0.49787 

0.50015 

0.49051 

0.48419 

0.48386 

0.48141 

0.47797 





0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 .00000 

0.00000 

0.00233 

0.00707 

0.01183 

0.01660 

0.02140 

0.02537 

0 . 02918 

0.03300 

0.03790 

0.04231 

0.04639 

0.04869 

0.05013 

0.04994 

0.05082 

0,05139 

0.05144 

0.05164 

0.05192 





1 0 

0 2 

2 2 

3 





STREAMLINE NUMBER 3 — 

STREAM FUNCTION « 0 

.0660 



1.40000 

1716.4800 

580.3772 

0.00307650 

.00061536 


0.000 

0 .00000 

301.1169 

26.5638 

0.120740 

0.0849164 





1.0000 

0.0010 







21 61 


81 20 

28 28 

1 1 

1 



0.0015377 

0.0038071 

52.3556 

1,7363 

13.0000 




0.000320 

0.002507 

0.012505 

0.022721 

0.053036 

0.043500 

0.054176 

0.064928 

0 . 075749 

0.086643 

0.097634 

0.108728 

0.116814 




0.0019246 

0.0076624 

0.0310414 

0.0494492 

0.0639724 

0.0756444 

0.0835645 

0.0877343 

0.0901731 

0.0912107 

0.0917784 

0.0921281 

0.0925331 




0.0015377 

0.0038071 

17.2296 

-1.8914 

12.0000 




0.001992 

0.012505 

0 . 022721 

0.033036 

0.043500 

0.054176 

0.064928 

0.075749 

0.086643 

0.097634 

0.108728 

0 .116803 





0.0030097 

0.0066846 

0.0200145 

0.0334604 

0. 0459736 

0.0574908 

0.0666530 

0 . 0731229 

0,0765852 

0.0780284 

0.0777994 

0.0773005 





-0.151661 

-0.119386 

-0.087776 

-0 . 056729 

-0.046853 

-0.036987 

-0.027134 

-0.017280 

-0.007403 

0.002507 

0.012505 

0.022721 

0.033036 

0.043500 

0.054176 

0 .064928 

0.075749 

0.086643 

0 .097634 

0.108728 

0.119859 

0.130821 

0.141635 

0 .162285 

0.182^88 

0.202545 

0.222573 

0.242572 





0.477011 

0.479328 

0.480943 

0.482355 

0.482990 

0.483824 

0.484708 

0.485999 

0.486423 

0.488575 

0.4900<«5 

0.490588 

0.492400 

0.493683 

0.494058 

0.494671 

0.495545 

0.497453 

0.498816 

0.499691 

0.499553 

0.500723 

0.500945 

0.501510 

0.501648 

0.501427 

0.501281 

0.501179 






0.0044123 0.0044123 0.0044123 0.0043910 0.0043778 0.0043678 0.0043544 0.0043436 

0.0043085 0.0042972 0.0042854 0.0042799 0.0042605 0.0042533 0.00^0468 0.0038859 

0,0037686 0.0037823 0.0037835 0.0037555 0.0037736 0.0036978 0.0037031 0.0037062 

0.0037222 0.0037203 0.0037305 0.0037453 

0.80894 0.78228 0.77002 0.76917 0.76916 0.76842 0.76795 0.76635 

0.76911 0.79603 0.S1327 0.83136 0.75130 0.68170 0.63448 0.60522 

0.59035 0.56646 0.55807 0.56667 0.56089 0.53216 0.53083 0.52955 

0.52646 0.52711 0.52545 0.52292 

0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 

0.00000 0.00104 0.00495 0.00913 0.01320 0.01750 0,02218 0.02592 

0.02957 0.03411 0.03841 0.04245 0.04553 0.04664 0.04653 0.04710 

0.04751 0.04751 0.04766 0.04788 

1 0 0 2 2 2 3 


X array 

RMSP 

BESP 


Printed Output List for Flow Solution 


QSONIC VERSION 1.5 SANZ SUPERCRITICAL STATOR HUB SECTION, FLOW SOLUTION 

CASE CONTROL INPUT ECHO 

CPARAMS 
NOFLOU= F 
RESTAR= F 
REMFSH= 4 

M3= 3, 7, 13, 25, 6^0 

NOZES= -9, -13, -25, -49, 6»0 

MESH1== 2 

MESHN= 4 

LAMDA0= 19.50 

CHORD’ 1.01460 

S= 0.8630 

NED= 88 

KN= 62 

PLE= 0.6999999E-02 
RTE= 0.290E-01 
THETL= 0,0 
THETT= 0 . 0 
CAMPER= 6 
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STABAC= 0.9990 
CHOP= 0.9899999 
SMOOTH= T 
LEOt<LY= F 
SLP1=: 1.0 
SLP2= 2.0 
SLP3= 1.0 
.SLP^»- 1.0 
MINF= 0.8010 
BETA1= 25.09999 
-BETA2= -15.60 
GAMMA= 1 .^0 

T0L5= 10000. 0» -0.1999999E-0<», 0 . 6 999999E-0^ , 0 . 5999999E-0<i , -0 . 999999<%E-06 
5^0.0 

OVEREL= 1.599999 
UNDERL= 1.0 
SUPPEL= 2.0 
MOWREL= 20 
NOTYET= 10 
TE€ARD= 2.0 
DAriP= 1.0 
CII= 0.^0 
IT= 30 
ALLOUT= T 
QUA.SI3= T 
NGTPM= 3 
RINF= 3.72^0 
WAKE= 0.0 
OriEGA= 0.0 
VAXIAL= 0.0 
FLOCO= 999.0 
M1NF2= 10.0 
BIGLAX= F 
F,END 


CALCULATED FLOW PARAHETERS 

ECONST 

CIP.CLN= 0.5^35562 
MINF2= 0.5237770 
B20B1= 0.8^88322 
R20R1= 1.050665 
RIHF= 3.72<«0 
RATLE= 3.810152 
CEMD 


RELAXATION BEGINS OH GRID OF 7 SURFACE CONTOURS INTERSECTED BY 25 RADIATING LINES 
SUPERSONIC RELAXATION PARArjETER. SUPREL= 0.2000DE 01 


ITERATION AVERAGE 
COUNT PHI CORRECTH 

30 0.l32^i39E-02 0 
60 0.303310E-03 0 
62 0.277992E-03 0 


RESIDUAL MAX MACH 
AT T. E. ON BLADE 

555933E-09 0.125951E 01 
295658E-09 0.129999E 01 
276595E-09 0.129591E 01 


CALCULATED RELATIVE CIRC 
CIRCULATION ERROR 

0.593598E 00 -0 . 196990E-09 
0.5935^5E 00 -0 .208348E-0<f 
0.593596E 00 -0 . 186^17E-0<i 


SUBSONIC SUPERSONIC 
RELAX. FAC. RELAX. FAC. 

0.160000E 01 0.200QOQE 01 
0 . 160000E 01 0 .200000E 01 
Q.160000E 01 0.200000E 01 


FINAL FLOW CALCULATION ON BLADE SURFACE FOR THE 7 BY 25 MESH 


PADIATING 



LINE NO. 

x/cx 

Y/CX 


Y 


STATIC 

MACH PRES COEF PHI 


XVEL 


YVEL DENSITY 


1 

2 
3 
9 

5 

6 
7 
S 
9 

10 

11 

12 

13 

19 

15 

16 

17 

18 

19 

20 
21 
22 
23 

25 


0.97785 

0.39962 

0.99928 

0.15395 

0.51052 

0.91197 

0.39329 

0.38739 

0.15221 

0.61619 

0.80899 

0.39353 

0.29187 

0.15299 

0.99291 

0.69070 

0.33126 

0.18160 

0.19100 

0.95173 

0.57196 

0.30299 

0.07092 

0.11969 

0.99687 

0.95978 

0.26005 

-0.03365 

0.07962 

0.95769 

0.35588 

0.20590 

-0.13051 

0.02919 

0.98869 

0.26005 

0.19617 

-0.21983 

-0.03153 

0.5^799 

0.17230 

0.08896 

-0.30169 

-0.08533 

0.61931 

0.09509 

0.09035 

-0.37360 

-0.13018 

0.67617 

0.03989 

0.00815 

-0.92977 

-0.16019 

0.61219 

0.00000 

0.00000 

-0.96229 

-0.16779 

0.93585 

0.00192 

0.02862 

-0.96096 

-0.19111 

0.72307 

0.09651 

0.08636 

-0.91889 

-0-08728 

1.11563 

0.12951 

0.16339 

-0.39618 

-0.01598 

1.29908 

0.22237 

0.29969 

-0.25996 

0.06026 

1.22505 

0.33096 

0.31337 

-0.15921 

0.12933 

1.09750 

0.99999 

0.35792 

-0.09795 

0.16538 

0.93931 

0.56086 

0.37939 

0.06057 

0.18582 

0.86319 

0.67965 

0.39208 

0.16669 

0.19769 

0.78239 

0.78219 

0.90003 

0.26689 

0.20511 

0.72767 

0.87899 

0.90979 

0.35660 

0.20999 

0.75593 

0.95866 

0.90728 

0.93138 

0.21187 

0.63819 

1.00000 

0.37980 

0.96992 

0.18158 

0.39361 

0.97785 

0.39962 

0.99928 

0.15395 

0.51055 


0.98029 

0.16129 

0.61095 

0.23723 

1.19972 

0.39359 

0.10687 

0.78732 

0.00838 

1.12330 

0.69899 

0.09208 

0.63689 

0.03060 

1.20099 

0.72513 

-0.02916 

0.57863 

0.09957 

1.22308 

0.73392 

-0.09013 

0.55991 

0.17399 

1.22568 

0.71976 

-0.15576 

0.59078 

0.29716 

1.22001 

0.65755 

-0.22207 

0.59772 

0.31797 

1.20301 

0.59281 

-0.29182 

0.59209 

0.38293 

1 . 16863 

0.39990 

-0.36526 

0.66567 

0.92809 

1.12369 

0.27331 

-0.93609 

0.79278 

0.92999 

1.08622 

0.39328 

-0.99230 

0.72736 

0.28786 

1.12320 

0.55523 

-0.51996 

0.36198 

0.92995 

1.17237 

0.06979 

-0.99851 

0.38572 

0.82056 

1.02227 

-0.67377 

-0.91775 

0.88261 

0.98592 

0.77309 

-0.99299 

-0.26777 

1.09891 

0.99877 

0.65606 

-0.87597 

-0.08865 

1.19507 

0.89663 

0.69990 

-0.59195 

0.07296 

1 .11577 

0.57819 

0.81936 

-0.30100 

0.20700 

1.09698 

0.31918 

0.90196 

-0.13959 

0.33012 

1 .05532 

0.16226 

0.95981 

0.03996 

0.93916 

0.97977 

0.09298 

1.01279 

0.16165 

0.53332 

0.91577 

0.05765 

1.05139 

0.09802 

0.61905 

0.99878 

0.03865 

1.03125 

0.25368 

0.68853 

0.78191 

-0.20995 

1.08012 

0.63995 

0.72139 

0.99973 

0.00196 

1 . 19761 

0.98027 

0.70978 

0.61099 

0.23736 

1.19971 


END OF CALCULATION ON GRID 1 OF 3 
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PELftXATIOH BEGINS ON GRID OF 13 SURFACE CONTOURS INTERSECTED BY 49 RADIATING LINES 


SUPERSONIC RELAXATION PARAMETER, SUPREL= 0.16000E 01 

ITERATION AVERAGE RESIDUAL MAX MACH CALCULATED RELATIVE CIRC SUBSONIC SUPERSONIC 

COUNT PHI CORRECTM AT T. E. ON BLADE CIRCULATION ERROR RELAX. FAC. RELAX. FAC. 

30 0.536023E-03 'Q . 250415E-04 0.111559E 01 0.543509E 00 -0 . 872868E-04 Q.1600Q0E 01 0.160000E 01 

60 0.141884E-03 0.774192E-05 0.113287E 01 0.543546E 00 -0 . 182030E-04 0.160000E 01 O.160Q00E 01 

90 0.946652E'04 0.729781E-05 0.114052E 01 0.543547E 00 -0 . 174354E-04 Q.160000E 01 0.160000E 01 

120 Q.727857E-04 0.505704E-05 0.114437E 01 0.343548E 00 -0 . 1 546 16E-04 0.160000E 01 0.160000E 01 

137 0.639313E-04 0.441199E-05 0.114672E 01 0.543549E 00 -0 . 14036 lE-04 0.160000E 01 0.160000E 01 


FINAL FLOW CALCULATION ON BLADE SURFACE FOR THE 13 BY 49 MESH 


RADIATING 
LINE NO. 

x/cx 

Y/CX 

X 

Y 

MACH 

STATIC 
PRES COEF 

PHI 

XVEL 

YVEL 

DENSITY 

1 

0.97795 

0.34310 

0.44928 

0.15345 

0.78982 

“0.04019 

0.15407 

0.92817 

0.32518 

0.98707 

2 

0.95042 

0.34133 

0.42350 

0.15180 

0.71379 

0.18381 

0.12828 

0.90083 

0.02841 

1.05830 

3 

0.91186 

0.34177 

0.38739 

0.15221 

0.58095 

0.47438 

0.09959 

0.74545 

-0.01310 

1.14793 

4 

0.86406 

0.34291 

0.34264 

0.15328 

0.51288 

0.61152 

0.06846 

0.66288 

-0.00317 

1.18926 

5 

0.80983 

0.34202 

0.29187 

0.15244 

0.47172 

0.68996 

0.03645 

0.61148 

0.02669 

1.21265 

6 

0.75174 

0.33791 

0.23748 

0.14860 

0.45030 

0.72904 

0.00405 

0.58220 

0.06097 

1.22424 

7 

0.69206 

0.32980 

0.18160 

0.14100 

0.44241 

0.74316 

-0.02851 

0.56712 

0.09797 

1.22841 

8 

0.63242 

0.31763 

0.12576 

0.12961 

0.44194 

0.74401 

-0.06128 

0.55881 

0.13522 

1.22866 

9 

0.57385 

0.30165 

0.07092 

0.11464 

0.44312 

0.74190 

-0.09405 

0.54991 

0.17274 

1.22804 

10 

0.51701 

0.28192 

0.01770 

0.09617 

0.44642 

0.73607 

-0.12663 

0.54117 

0.21014 

1.22631 

11 

0.46216 

0.25890 

-0.03365 

0.07462 

0.45311 

0.72415 

“0.15922 

0.53505 

0.24601 

1.22279 

12 

0.40938 

0.23310 

-0.08308 

0.05046 

0.46287 

0.70651 

-0.19199 

0.53152 

0.28066 

1.21756 

13 

0.35872 

0.20499 

-0 .13051 

0.02414 

0.47891 

0.67697 

“0.22507 

0.53551 

0.31444 

1.20879 

14 

0.31005 

0 .17550 

-0.17608 

-0.00347 

0.50525 

0.62701 

“0.25894 

0.55286 

0.34850 

1.19390 

15 

0.26332 

0.14553 

-0.21983 

-0.03153 

0.53928 

0.56047 

“C .29403 

0.58172 

0.38067 

1.17394 

16 

0.21849 

0.11613 

-0.26180 

-0.05906 

0.57585 

0.48657 

-0.33010 

0.61656 

0.40830 

1.15163 

17 

0.17595 

0.08807 

-0.30164 

-0.08533 

0.61045 

0.41480 

“0.36652 

0.65401 

0.42679 

1.12978 

18 

0 .13592 

0.06249 

-0.33911 

-0 ,10928 

0.63735 

0.35791 

“0.40215 

0.68928 

0.43079 

1.11235 

19 

0.09908 

0.04017 

-0 . 37360 

-0.13018 

0.61765 

0.39878 

“0.43549 

0.67922 

0.40239 

1.12488 

20 

0.06650 

0.02150 

-0.40411 

-0.14766 

0.59873 

0.43880 

-0.46200 

0.67744 

0.55953 

1.13710 

21 

0.03910 

0 .00812 

-0.42977 

-0.16019 

0.17308 

-0.01450 

“0.48418 

-0.21084 

“0.04639 

0.99534 

22 

0.01825 

0 .00003 

-0.44929 

-0 . 16776 

6.22106 

-1.54228 

-0.48276 

6.55885 

-0.20155 

0.43053 

23 

0.00441 

0 .00000 

-0.46224 

-0 . 16779 

56 . 36440 

-2.22657 

-0.62053 

4.47947 

5.98375 

0.00001 

24 

0.00000 

0 .00997 

-0.46638 

-0.15846 

6 . 32647 

-1 .36754 

-0.49262 

-2.48300 

6.43079 

0.50647 

25 

0.00632 

0.02849 

-0.46046 

-0 . 14111 

0.80365 

-0.02053 

-0.47991 

0.44348 

0.89849 

0.99340 

26 

0.02331 

0.05436 

-0 . 44455 

-0.11689 

0.95377 

“0 . 34364 

-0.44982 

0.70725 

0.92268 

0.88715 

27 

0 .05072 

0.08598 

-0.41889 

-0 . 08728 

1 . 08899 

-0.61440 

-0.39923 

0.87627 

0.95791 

0.79403 

28 

0 . 08630 

0 .12269 

“0 . 38557 

-0.05291 

1 .1 0408 

-0 .64386 

-0.33654 

0.93137 

0.92517 

0.78364 

29 

0.12638 

0 .16267 

-0 . 34618 

-0.01548 

1.11713 

-0.66942 

-0.26475 

0.97728 

0.89504 

0.77458 

30 

0.17532 

0.20375 

-0 . 30223 

0 . 02298 

1 .13097 

-0.69652 

-0 . 18676 

1.02708 

0.85795 

0.76493 

31 

0.22580 

0.24357 

-0.25496 

0.06026 

1.14337 

-0.72066 

-0.10561 

1 . 08273 

0.80612 

0.75629 

32 

0.27873 

0.28022 

-0.20540 

0.09458 

1 .14589 

-0.72614 

-0.02366 

1 .13823 

0.72986 

0.75432 

33 

0.33341 

0.31199 

-0.15421 

0 .1 2433 

1.11188 

-0.66160 

0.05650 

1 .16989 

0.61119 

0.77736 

34 

0.38953 

0.33744 

-0.10166 

0.14815 

1.04958 

-0.53805 

0.13123 

1 .17278 

0.45944 

0.82071 

35 

0.44689 

0.35584 

-0 . 04795 

0 .16538 

0.98921 

-0.41184 

0.20072 

1.15718 

0.31549 

0.86408 

36 

0.50494 

0.36841 

0.00639 

0 .17716 

0.92702 

-0.27760 

0.26586 

1.11526 

0.21356 

0.90927 

37 

0.56280 

0.37766 

0.06057 

0.18582 

0.87224 

-0.15764 

0.32627 

1.06651 

0.15402 

0.94890 

38 

0.61998 

0 . 38477 

0.11411 

0 .19247 

0.82790 

-0.05995 

0.38319 

1.02295 

0.11649 

0.98070 

39 

0.67608 

0 .39035 

0.16664 

0 .19769 

0.78069 

0.04446 

0.43622 

0.97327 

0.08971 

1 .01422 

'•O 

ft . 7 3 0 r, 

ft . 39ft78 

0.2177? 

0 . POl 84 

0 . 719)9 

0 r 136 OB 

D .4A493 

0 .92809 

0 .07000 

1 . 04310 

41 

0.78310 

0.39827 

0.26684 

0.20511 

0.71222 

0.19585 

0.52978 

0.89821 

0.05561 

1.06207 

42 

0.83280 

0.40099 

0 .31337 

0.20765 

0.69987 

0.22296 

0.57121 

0.88464 

0 . 04403 

1.07054 

43 

0.87897 

0.40295 

0.35660 

0.20949 

0.70178 

0.21876 

0.60940 

0.88730 

0.03367 

1.06923 

44 

0.92086 

0.40428 

0.39582 

0.21073 

0.75784 

0.09340 

0.64445 

0.95119 

0.03089 

1.02979 

45 

0.95884 

0.40549 

0.43138 

0.21187 

0.80160 

-0.02587 

0.68025 

0.99022 

-0.13275 

0.99169 

46 

0.99107 

0.39581 

0.46156 

0.20281 

0.53391 

0.40780 

0.71232 

0.57170 

-0.37326 

1.12764 

47 

1.00000 

0.37315 

0.46992 

0.18158 

0.12115 

0.92023 

0.72139 

0.15326 

“0.04199 

1.28029 

48 

0 . 99231 

0.35274 

0.46272 

0.16248 

0.53366 

0.43922 

0.71445 

0.49221 

0.47437 

1.13723 

49 

0.97795 

0.34310 

0.44928 

0.15345 

0.78988 

-0.04034 

0.69762 

0.92820 

0.32530 

0.98702 


END OF CALCULATION ON GRID 2 OF 3 


RELAXATION BEGINS ON GRID OF 25 SURFACE CONTOURS INTERSECTED BY 97 RADIATING LINES 

SUPERSONIC RELAXATION PARAMETER, SUPREL= 0.12000E 01 

ITERATION AVERAGE RESIDUAL MAX MACH CALCULATED RELATIVE CIRC SUBSONIC SUPERSONIC 

COUNT PHI CORRECTN AT T. E. ON BLADE CIRCULATION ERROR RELAX. FAC. RELAX. FAC. 

30 0.158286E-03 -0 , 83157 lE-05 0.113079E 01 0.543551E 00 -0 . 1 00884E-04 0.160000E 01 0.120000E 01 

60 0.635203E-04 0.650720E-06 0.115129E 01 0.543554E 00 -0 . 405730E-05 0.160000E 01 0.120000E 01 

64 0.589319E-04 0.706636E-06 0.115269E 01 0.543555E 00 -0 . 219314E-05 0.160000E 01 0.120000E 01 


FINAL FLOW CALCULATION OH BLADE SURFACE FOR THE 25 BY 97 MESH 


RADIATING 
LINE NO. 

X/CX 

Y/CX 

X 

Y 

MACH 

STATIC 
PRES COEF 

PHI 

XVEL 

YVEL 

DENSITY 

1 

0.97795 

0.34420 

0.44928 

0.15345 

0.96863 

-0.38816 

0.15264 

1.12082 

0.35810 

0.87212 

2 

0.96587 

0.34232 

0.43796 

0.15169 

0.82943 

“0.07382 

0.13970 

1.02698 

0.08536 

0.97621 

3 

0.95043 

0.34244 

0.42350 

0.15180 

0.68838 

0.24513 

0.12618 

0.87228 

-0.00775 

1.07745 

4 

0.93244 

0.34262 

0.40666 

0.15197 

0.61172 

0.41190 

0.11257 

0.78239 

“0.00897 

1.12889 

5 

0.91186 

0.34287 

0.38739 

0.15221 

0.57356 

0.49132 

0.09799 

0.73662 

-0.01340 

1.15307 

6 

0.88898 

0.34341 

0.36596 

0.15271 

0.53887 

0.56129 

0.08263 

0.69450 

“0.01695 

1.17419 

7 

0.86407 

0.34401 

0.34264 

0.15328 

0.50585 

0.62590 

0.06694 

0.65423 

“0.00805 

1.19357 

8 

0.83758 

0.34399 

0.31783 

0.15326 

0.48168 

0.67189 

0.05116 

0.62435 

0.01014 

1.20728 

9 

0.80985 

0.34312 

0.29187 

0.15244 

0.46645 

0.70018 

0.03523 

0.60493 

0.02688 

1.21568 

10 

0.78114 

0.34152 

0.26498 

0.15094 

0.45516 

0.72069 

0.01913 

0.58995 

0.04232 

1.22176 

11 

0.75176 

0.33902 

0.23748 

0.14860 

0.44655 

0.73612 

0.00298 

0.57761 

0.05994 

1.22633 

12 

0.72200 

0,33549 

Q .20961 

0.14529 

0.44149 

0,74512 

“0.01320 

0.56895 

0.07875 

1.22899 

13 

0.69209 

0.33090 

0.18160 

0.14100 

0.45979 

0.74812 

“0.02943 

0.56386 

0.09766 

1.22987 

14 

0.66219 

0.32531 

0.15360 

0.13576 

0.44061 

0.74669 

-0.04574 

0.56132 

0.11647 

1.22945 

15 

0.63245 

0.31874 

0.12576 

0.12961 

0.44223 

0,74381 

“0.06214 

0.55916 

0.13534 

1.22860 

16 

0.60299 

0.31121 

0.09817 

0.12256 

0.44269 

0.74301 

“0.07856 

0.55499 

0.15370 

1.22836 

17 

0.57389 

0.30275 

0.07092 

0.11464 

0.44233 

0.74363 

“0.09491 

0.54911 

0.17205 

1.22855 


45 



18 

0.59521 

0.29336 

0.09907 

0.10585 

0.99329 

0.79201 

-0.11119 

0.59385 

0.19197 

1.22807 

19 

0.51705 

0.28302 

0.01770 

0.09617 

0.99585 

0.73790 

-0.12795 

0.59020 

0.21069 

1.22671 

20 

0.98938 

0.27188 

-0.00821 

0.08573 

0.99919 

0.73195 

-0.19373 

0.53799 

0.22859 

1.22995 

21 

0.96221 

0.26001 

-0.03365 

0.07962 

0.95313 

0.72938 

-0.16005 

0.53500 

0.29619 

1.22286 

22 

0.93556 

0.29793 

-0.05860 

0.06289 

0.95765 

0.71623 

-0.17691 

0.53289 

0.26372 

1,22099 

23 

0.90992 

0.23921 

-0.08308 

0.05096 

0.96303 

0.70699 

-0.19289 

0.53126 

0.28160 

1.21759 

29 

0.58387 

0.22036 

-0.10701 

0.03799 

0.96991 

0.69980 

-0.20935 

0.53120 

0.29829 

1.21909 

25 

0.35877 

0.20610 

-0.13051 

0.02919 

0.97723 

0.68035 

-0.22598 

0.53353 

0.31371 

1.20979 

26 

0.33919 

0.19199 

-0.15353 

0.01096 

0.98892 

0,65892 

-0.29275 

0.59010 

0.33091 

1.20327 

27 

0.31010 

0.17661 

-0.17608 

-0.00397 

0.50915 

0.62997 

-0.25983 

0.55139 

0.39831 

1.19963 

28 

0.28698 

0.16166 

-0.19820 

-0.01797 

0.52191 

0.59602 

-0.27721 

0.56519 

0.36606 

1.18962 

29 

0.26338 

0.19669 

-0.21983 

-0.03153 

0.59007 

0.55922 

-0.29992 

0.58166 

0.38259 

1.17357 

30 

0.29068 

0.13180 

-0.29108 

-0.09593 

0.55823 

0.52281 

-0.31299 

0.59918 

0.39603 

1.16259 

31 

0.21856 

0.11725 

-0.26180 

-0.05906 

0.57589 

0.98692 

-0.33105 

0.61621 

0.90883 

1,15173 

32 

0.19699 

0.10301 

-0.28200 

-0.07239 

0.59939 

0.99869 

-0.39929 

0.63975 

0.^2108 

1.14012 

33 

0.17601 

0.08919 

-0.30169 

-0.08533 

0.61222 

0.91132 

-0.36795 

0.65500 

0.92912 

1.12872 

39 

0.15565 

0.07606 

-0.32071 

-0.09763 

0.62731 

0.37993 

-0.38597 

0.67398 

0.93278 

1.11895 

35 

0.13599 

0.06361 

-0.33911 

-0.10928 

0.63971 

0.35300 

-0.90312 

0.69082 

0.93357 

1.11089 

36 

0.11713 

0.05198 

-0.35678 

-0,12018 

0,69771 

0.33590 

-0.92020 

0.70992 

0.92867 

1.10557 

37 

0.09916 

0.09129 

-0.37360 

-0.13018 

0.65213 

0.32696 

-0.93693 

0.71997 

0.92289 

1.10266 

38 

0.08226 

0.03192 

-0,38993 

-0.13993 

0.65852 

0.31265 

-0.95170 

0.72635 

0.91732 

1.09890 

39 

0.06658 

0.02262 

-0.90911 

-0.19766 

0.65601 

0.31770 

-0.96590 

0,73689 

0.39215 

1.09996 

90 

0.05209 

0.01535 

-0.91773 

-0.15998 

0.69073 

0.35079 

-0.97850 

0.73385 

0.35869 

1.11019 

91 

0.05918 

0.00929 

-0.92977 

-0.16019 

0.61669 

0.90173 

-0.98926 

0.71530 

0.33129 

1,12579 

9 2 

0.02795 

0.00920 

-0.99028 

-0.16992 

0.60532 

0.92370 

-0.99813 

0.72395 

0.27592 

1.13250 

93 

0.01833 

0.00116 

-0.99929 

-0.16776 

0.51533 

0.59251 

-0.50552 

0.69629 

0.15860 

1.18357 

99 

0.01053 

0.00000 

-0.95659 

-0.16885 

0.90030 

0.80589 

-0.50959 

0.52226 

-0.00576 

1.29687 

95 

0.00999 

0.00115 

“0.96229 

-0 . 16779 

0.22327 

0.99709 

-0.51299 

0.28446 

-0.07376 

1.30255 

96 

0.00000 

0.00957 

-0.96695 

-0.16957 

0.63265 

-0.19069 

-0.51296 

0.09558 

0.78117 

0.95995 

97 

0.00008 

0.01109 

-0.96638 

-0.15896 

0.99629 

-0.90129 

-0.50322 

0 . 09190 

1.19967 

0.86768 

98 

0.00209 

0.01938 

-0.96950 

-0 .15070 

0.69529 

0.20361 

-0.99809 

0.28935 

0.82988 

1.06450 

99 

0.00690 

0.02962 

-0.96096 

-0.19111 

0.99703 

-0.32675 

-0.98689 

0.52989 

1.02978 

0.89283 

50 

0 .01355 

0.09171 

-0.95376 

-0,12979 

1.07979 

-0.58790 

-0.97050 

0.70168 

1.07568 

0.80351 

51 

0. 02339 

0. 05598 

-0.99955 

-0.11689 

1.12790 

-0.69157 

-0.99998 

0.81752 

1.05573 

0.76670 

52 

0.03598 

0.07098 

-0.93276 

-0.10285 

1.13691 

-0,70652 

-0.92976 

0.87787 

1.01695 

0.76135 

53 

0.05079 

0.08710 

-0.91889 

-0.08728 

1.12769 

-0.68966 

-0.39678 

0.90269 

0.98380 

0.76738 

59 

0.06766 

0.10986 

-0.90310 

-0.07065 

1.11558 

-0.66606 

-0.36622 

0.91955 

0.95228 

0.77578 

55 

0.08658 

0.12381 

-0.38557 

-0.05291 

1.10170 

-0.63886 

-0.33333 

0.92922 

0.92916 

0.78591 

56 

0.10669 

0.19351 

-0.36655 

-0.03997 

1.08801 

-0.61179 

-0.29870 

0.93893 

0,89587 

0,79997 

57 

0.12895 

0.16378 

-0.39618 

-0.01598 

1.08058 

-0.59690 

-0.26269 

0.95013 

0.87287 

0.80018 

58 

0.15137 

0.18936 

-0.32971 

0.00379 

1.08532 

-0.60693 

-0.22590 

0.97219 

0.85518 

0.79689 

59 

0.17538 

0.20986 

-0.30223 

0.02298 

1.09602 

-0.62768 

-0.18689 

1.00111 

0.83728 

0.78936 

60 

0.20027 

0.22509 

-0.27893 

0.09187 

1 .10699 

-0.69932 

-0.19799 

1.03069 

0.81765 

0.78171 

61 

0.22586 

0.29967 

-0.25996 

0 . 06026 

1 . 12032 

-0-67593 

-0.10750 

1.06332 

0.79601 

0.77299 

62 

0.25206 

0.26355 

-0.25093 

0.07793 

1.13612 

-0.70611 

-0.06701 

1.10190 

0.76879 

0.76150 

63 

0.27879 

0.28133 

-0.20590 

0.09958 

1.19889 

-0.73059 

-0.02629 

1.13983 

0.73287 

0.75272 

69 

0.30599 

0.29789 

-0.37998 

0.11009 

1.15275 

-0.73801 

0.01931 

1.17189 

0.68753 

0.75006 

65 

0.33396 

0.31309 

-0.15921 

0.12933 

1 .13592 

-0.70506 

0.05990 

1.18752 

0.62610 

0.76188 

66 

0.36132 

0.32677 

-0.12812 

0.13713 

1.08696 

-0.60979 

0.09306 

1.17751 

0,59082 

0.79567 

67 

0.38958 

0.33859 

-0.10166 

0.19815 

1.03567 

-0.50682 

0.12939 

1.16209 

0.95002 

0.83153 

68 

0.91816 

0.39850 

-0.07990 

0.15798 

1.02099 

-0.97509 

0.16929 

1.17095 

0.37988 

0.89298 

69 

0.99699 

0.35699 

-0.09795 

0 .16538 

1.00829 

-0.99991 

0.19895 

1.17768 

0.31386 

0.85109 

70 

0.97596 

0.36370 

-0.02078 

0.17171 

0.96916 

-0.35661 

0.23255 

1.19668 

0. 25205 

0.88278 

71 

0.50998 

0.36951 

0.00639 

0.17716 

0.92239 

-0.26669 

0.26929 

1.11112 

0.20986 

0-91292 

72 

0.53395 

0.37950 

0.03352 

0 . 18183 

0.89571 

-0-20867 

0.29505 

1.08891 

0.17699 

0.93212 

73 

0.56283 

0.37876 

0.06057 

0.18582 

0.87026 

-0.15290 

0.32981 

1.06977 

0.15199 

0.95095 

79 

0.59152 

0.38258 

0.08793 

0 . 18939 

0.85029 

-0.10888 

0.35362 

1.09562 

0.13189 

0.96982 

75 

0.62001 

0.38587 

0.11911 

0.19297 

0.83132 

-0.06712 

0.38169 

1.02689 

0.11512 

0.97837 

76 

0.69822 

0.38883 

0,19053 

0.19525 

0.80632 

-0.01193 

0.90883 

1.00076 

0.10121 

0.99617 

77 

0.67611 

0.39195 

0.16669 

0.19769 

0.77829 

0.05010 

0.93980 

0.97071 

0.08857 

1.01602 

78 

0.70361 

0.39380 

0,19239 

0.19990 

0.75379 

0.10932 

0.95960 

0.99901 

0.07819 

1.03325 

79 

0.73066 

0.39588 

0.21772 

0.20189 

0.73955 

0.19681 

0.98337 

0.92286 

0.06909 

1.09666 

80 

0.75720 

0.39773 

0.29257 

0.20358 

0.71929 

0.18090 

0.50618 

0.90596 

0.06191 

1.05722 

81 

0.78312 

0.39936 

0.26689 

0.20511 

0.70721 

0.20691 

0.52807 

0.89250 

0.05992 

1.06553 

82 

0.80836 

0.90081 

0.29097 

0.20696 

0.69873 

0.22596 

0.59910 

0.88305 

0,09933 

1.07133 

S3 

0.83281 

0.90208 

0.31337 

0.20765 

0.69955 

0.23959 

0.56929 

0.87851 

0.09399 

1.07917 

89 

0.85639 

0.90317 

0.33595 

0.20867 

0.69350 

0.23689 

0.58871 

0.87756 

0.03820 

1 . 07989 

85 

0.87898 

0.90905 

0.35660 

0.20999 

0.69929 

0.23527 

0.60731 

0.87866 

0.03228 

1.07938 

86 

0.90096 

0.90975 

0.37672 

0.21015 

0.69751 

0.22819 

0.62502 

0.88256 

0.02853 

1 . 07216 

87 

0.92086 

0.90537 

0.39582 

0.21073 

0.70996 

0.20192 

0.69195 

0.89631 

0.02866 

1.06397 

88 

0.99023 

0.90599 

0.91396 

0.21131 

0.75597 

0.09853 

0.65892 

0.99919 

0.03092 

1.03191 

89 

0.95889 

0.90658 

0.93138 

0.21187 

0.90915 

-0.23525 

0.67570 

1.11060 

-0.02671 

0.92339 

90 

0.97679 

0.90537 

0.99819 

0.21073 

0.91662 

-0.31289 

0.69630 

1.07323 

-0.31976 

0.89798 

91 

0.99107 

0.39691 

0.96156 

0.20281 

0.65316 

0.25183 

0.71300 

0.62951 

-0.59368 

1.07959 

92 

0 . 99791 

0-38618 

0.96797 

0 .19276 

0.30668 

0.89273 

0.72050 

0.18193 

-0.35822 

1.27229 

93 

1.00000 

0.37929 

0.96992 

0.18158 

0.09073 

1.10179 

0.72230 

0.03883 

0.03732 

1.33261 

99 

0.99799 

0.36307 

0,96799 

0.17112 

0.31969 

0.90129 

0.71992 

0.17359 

0.38127 

1.27978 

95 

0.99231 

0.35389 

0.96272 

0.16298 

0.57935 

0-99651 

0.71398 

0.98180 

0.56985 

1,13995 

96 

0 . 98586 

0.39790 

0.95668 

0.15691 

0.86911 

-0.17913 

0.70692 

0.89370 

0.58167 

0.99399 

97 

0.97795 

0.39920 

0.94928 

0.15345 

0. 96866 

-0.38822 

0.69619 

1,12089 

0.35812 

0.87210 


END OF CALCULATION ON GRID 3 OF 3 


CALCULATED FLOW FIELD AT ALL GRID POINTS 

SURFACE CONTOUR NO. 1 = PERIODIC BOUNDARY 

SURFACE CONTOUR HO. 25 = BLADE SURFACE 


RADIATING LINE NO. 1 BEGINS AT DOWNSTREAM INFINITY 
AND IS REPEATED AS LINE NO. 97 




SURFACE 








FLOW 

STATIC 



CONTOUR 

NUMBER 

X/CX 

Y/CX 

PHI 

XVEL 

YVEL 

DENSITY 

MACH 

ANGLE 

PRES COEF 

RADIATING 

LINE 











NUMBER 

1 

1 

2.61888 

0.28898 

1.98505 

-0.29396 

0.93997 

0.51703 

0.52378 

119.31993 

-1.55813 



2 

1.79953 

0.28823 

0.93527 

1.39087 

-0.27375 

0.61525 

0.52378 

-11.53872 

3.77358 



3 

1.59681 

0.28861 

0.53292 

0.68003 

0.09605 

1.18072 

0.52378 

3.87389 

0.62970 



9 

1.48078 

0.28922 

0.95907 

0.68019 

0.09635 

1.18063 

0.52825 

3.89807 

0.58273 



5 

1.39876 

0.28985 

0.90689 

0.67979 

0.09673 

1.18089 

0.52791 

3.93293 

0.58591 



6 

1.33567 

0.29069 

0.36680 

0.67831 

0.09739 

1.18150 

0.52678 

3.99231 

0.58563 



7 

1.28975 

0.29159 

0.33955 

0.67569 

0.09828 

1.18273 

0.52969 

9.08721 

0.58971 



8 

1.29228 

0.29272 

0.30780 

0.67219 

0.09969 

1.18938 

0.52189 

9.22799 

0.59520 



9 

1.20611 

0.29909 

0.28516 

0.66779 

0.05172 

1.18637 

0.51898 

9.92880 

0.60185 



10 

1.17982 

0.29559 

0.26579 

0.66279 

0.05962 

1.18865 

0.51956 

9,71101 

0.60996 



11 

1.19797 

0.29721 

0.29892 

0.65708 

0.05866 

1.19119 

0.51025 

5.10199 

0.61779 


46 





12 

1.12338 

0.29909 

0.23627 

0.65087 

0,06622 

1.19378 

0.50566 

5.63677 

0.62663 



13 

1.10202 

0.30116 

0.22166 

0.66613 

0.07179 

1.19652 

0.50085 

6.35969 

0.63580 



l<t 

1.08306 

0.30366 

0.21026 

0.63660 

0,08173 

1.19961 

0.69576 

7.31591 

0.66568 



15 

1.06616 

0.30592 

0.20063 

0 . 62887 

0.09662 

1.20208 

0.69098 

8.55638 

0.65666 



16 

1.05111 

0.30862 

0.19187 

0.62Z96 

O.UlU 

1.20369 

0.68863 

10.11721 

0.65917 



17 

1.03776 

0.31153 

0.18663 

0.62088 

0.13178 

1.20257 

0.68999 

11.98352 

0.65608 



18 

1.02588 

0.31668 

D. 17796 

0.62670 

0.15686 

1.19820 

0.69759 

16.09326 

0.66163 



19 

1.01563 

0 .31806 

0.17235 

0.63717 

0.18630 

1.18885 

0.51366 

16.29860 

0.61013 



20 

1 .00629 

0.32170 

0.16752 

0.66185 

0.21956 

1.17267 

0.56105 

18,35106 

0.55558 



21 

0.99839 

0.32559 

0.16338 

0.70310 

0.25508 

1.16661 

0.58295 

19.96038 

0.66939 



22 

0.99166 

0.32976 

0.15985 

0,76585 

0.29026 

1.10736 

0.66278 

20.75589 

0,36162 



23 

0.98601 

0.33621 

0.15689 

0.85651 

0.32106 

1.05156 

0.72386 

20.59166 

0.16237 



2<4 

0.98166 

0.33896 

0.15665 

0 . 97207 

0.36396 

0.97523 

0.83009 

19.68599 

-0.07682 



25 

0.97795 

0.36620 

0.15266 

1.12082 

0.35810 

0.87212 

0.96863 

17.71860 

-0.38816 

RADIATING 

LINE 











NUMBER 

2 













1 

1.66562 

-0.16596 

0.53726 

0.12880 

1.56907 

0.86221 

0.52378 

85.26683 

8.36990 



2 

1.59936 

-0.00776 

0.61053 

1.01730 

0.15876 

0.91636 

0.52378 

8.86996 

2.83198 



3 

1.51637 

0.09003 

0.67167 

0.68116 

0 . 06697 

1 . 18015 

0.52378 

3.96507 

0.63811 



6 

1 .63676 

0.16506 

0.62317 

0.68186 

0.06860 

1.17975 

0.52975 

6.07679 

0.57978 



5 

1.36799 

0.17852 

0.38213 

0.68169 

0.06968 

1.17989 

0.52951 

6.16959 

Q. 58026 



6 

1.31253 

0.20088 

0.36783 

0.68026 

0 .05106 

1.18063 

0.52859 

6.29069 

0.58206 



7 

1.26578 

•0.21705 

0.31890 

0.67817 

0.05283 

1.18137 

0.52700 

6.65635 

0,58520 



8 

1.22582 

0.22969 

0.29622 

0.67565 

0.05518 

1.18259 

0.52693 

6,67076 

0.58926 



9 

1.19126 

0.23957 

0.27295 

0.67228 

0.05822 

1.18398 

0.52256 

6.96971 

0.59388 



10 

1 .16101 

0.26813 

0.25666 

0.66886 

0.06216 

1.18565 

0.52005 

5.30820 

0.59877 



11 

1.13637 

0.25563 

0.23828 

0.66530 

0.06717 

1.18689 

0.51757 

5.76685 

0.60359 



12 

1.11076 

0.26262 

0.22607 

0.66176 

0.07362 

1 . 18826 

0.51526 

6.36819 

0.60811 



13 

1.08976 

0.26876 

0.21155 

0.65797 

0.03168 

1 . 18957 

0.51295 

7.07636 

0.61253 



1<» 

1 .07099 

0.27677 

0.20052 

0.55620 

0.09166 

1 .19072 

0.51096 

7.95851 

0.61637 



15 

1.05625 

0.28059 

0.19083 

D. 65186 

0.10337 

1.19097 

0.51051 

9.01066 

0.61721 



16 

1.03931 

0.28631 

0.18231 

0.65289 

0.11752 

1.18936 

0.51325 

10.20358 

0.61183 



17 

1.02598 

0.29200 

0.17682 

0.65892 

0.13361 

1.18507 

0.52056 

11.66280 

0.59752 



18 

1.01616 

0.29773 

0.16826 

0.67176 

0.15105 

1.17721 

0.53379 

12.67282 

0.57132 



19 

1.00371 

0.30353 

0.16265 

0.69365 

0. 16863 

1 .16667 

0.55667 

13.66786 

0.52970 



20 

0.99655 

0 . 30966 

0.15736 

0.72602 

0 . 18626 

1 . 16660 

0.58380 

16.23916 

0.66936 



21 

0.98660 

0.31555 

0.15287 

0.77070 

0 . 19666 

1 . 12161 

0.62227 

16.17536 

0.38809 



22 

0 . 97981 

0 . 32185 

0.16891 

0 . 82697 

0. 1 9567 

1 .09039 

0.66901 

13.31188 

0,28675 



23 

0.97612 

0.32838 

0 . 16562 

0.89162 

0.18296 

1 .05616 

0.72162 

11.59687 

0.17063 




0.96969 

0.3351 9 

0.16235 

0 . 95869 

0 .15661 

1.01579 

0.77525 

9.16329 

0.06937 



25 

0.96587 

0.56232 

0 .13970 

1 .02698 

0 . 08536 

0.97621 

0.82963 

6.75160 

-0.07382 

RADIATING 

LINE 











HUMBER 

3 













1 

1.66305 

-0.16667 

0.61391 

0.68296 

0 . 06555 

1 .17932 

0.52378 

3.81530 

0.65056 



2 

1.62930 

-0.07861 

0 .60905 

0 .68316 

0 . 06679 

1.17918 

0.52378 

3.91866 

0.65258 



3 

1.39537 

-0.00360 

0 . 39092 

0 .68666 

0 . 06876 

1 .17869 

0.52378 

6.07682 

0.66280 




1 . 35250 

0 .0566 0 

0 . 36620 

0 .68520 

0 . 05066 

1.17806 

0 .53260 

6.22870 

0 . 57615 



5 

1 .30852 

0 .09768 

0.36011 

0 .68679 

0 .052 30 

1 .17820 

0.53236 

6.36772 

0.57662 



6 

1 .267 06 

0.13005 

0.31520 

0 . 68376 

0.05635 

1.17863 

0.53166 

6.56691 

0.57606 



7 

1.22922 

0.15525 

0,29236 

0.68220 

0 . 05686 

1 .17927 

0.53055 

6.76661 

0.57820 



8 

1 .19517 

0.17561 

0.27177 

0 . 68035 

0.05991 

1 . 18006 

0.52925 

5.03197 

0.58076 



9 

1 .16663 

0 .19206 

0.25333 

0.67839 

0 .06358 

1 . 18082 

0.52792 

5.35656 

0.58336 



10 

1 .1 3726 

0.20618 

0.23685 

0 .6 7667 

0.06802 

1.18156 

0.52671 

5.76230 

0.58573 



11 

1.11265 

0.21852 

0 , 22213 

0.67665 

0 .07360 

1 . 18213 

0.52569 

6.20898 

0.58773 



12 

1 .09057 

0.22956 

0.20900 

0.67269 

0 .07976 

1.18272 

0 . 52669 

6.76016 

0.58970 



13 

1.07071 

0.23960 

0.19730 

0.67067 

0.08697 

1.18326 

0.52377 

7.38871 

0.59169 



U 

1 .05285 

0 . 26898 

0.18691 

0.66969 

0 .09526 

1.18319 

0.52389 

8.09608 

0.59126 



15 

1 .03678 

0.25785 

0.17767 

0 .67126 

0.10651 

1.18178 

0.52628 

8.86960 

0.58653 



16 

1.02235 

0.26637 

0.16966 

0 .67663 

0 .11625 

1.17865 

0.53189 

9.58398 

0.57565 



17 

1 .00966 

0.27667 

0.16216 

0.68687 

0.12373 

1 .17269 

0.56151 

10.21190 

0.55629 



IS 

0.99793 

0.28285 

0 . 1 5566 

0.70296 

0.13183 

1 .16606 

0-55575 

10.62161 

0.52761 



19 

0.98773 

0.29098 

0.16987 

0 .72539 

0 .13687 

1 .15228 

0.57676 

10.68520 

0.68872 



20 

0 . 97876 

0.29916 

0 . 16669 

0.75360 

0.13663 

1.13768 

0.59785 

10.27616 

0.66070 



21 

0 . 97096 

0.30760 

0 .16005 

0.78567 

0.12868 

1 .12123 

0.62327 

9.30166 

0.38687 



22 

0 . 96627 

0.31580 

0.13591 

0.81808 

0.11111 

1 .10671 

0.66825 

7.73656 

0.3331 1 



23 

0 . 95863 

0.32660 

0.13223 

0.86609 

0.08353 

1 .09051 

0.66931 

5.63853 

0 .28716 



2<i 

0 .95605 

0.33322 

0.12901 

0.86699 

0 .06763 

1 .08098 

0.68318 

3.13862 

0.25666 



25 

0.95063 

0.56266 

0.12618 

0.87228 

-0.00775 

1.07765 

0.68838 

-0.50906 

0.26513 

RADIATING 

LINE 











NUMBER 


1 

1 . 32367 

-0.16789 

0.33715 

0.68865 

0.06696 

1.17657 

0.53509 

3.73567 

0.56921 



2 

1 .31638 

-0.10591 

0.33533 

0.68897 

0.06753 

1.17633 

0.53550 

3.96673 

0 . 56861 



3 

1.29861 

-0.06851 

0.32653 

0.68975 

0.05005 

1 .17587 

0.53628 

6.15022 

0.56687 




1 .27356 

0.00166 

0.31278 

0.68989 

0.05213 

1.17572 

0.53652 

6.32100 

0 . 56638 



5 

1 .26679 

D .06560 

0.29635 

0.68950 

0 . 05656 

1.17583 

0.53635 

6.52267 

0.56672 



6 

1 .21505 

0.07816 

0 .27902 

0.68873 

0.05725 

1.17609 

0.53591 

6.75160 

0.56760 



7 

1 . 18592 

0 .10706 

0.26188 

0.68776 

0.06029 

1 .17666 

0.53531 

5.01013 

0,56877 



8 

1 . 15829 

0.13131 

0.26556 

0.68668 

0.06375 

1.17681 

0.53670 

5.30383 

0.56999 



9 

1 .13269 

0.15201 

0.23027 

0.68560 

0 , 06770 

1 . 17715 

0.53613 

5.63968 

0.57112 



10 

1 .10866 

0 . 17001 

0.21618 

0.63666 

0.07216 

1.17769 

0.53355 

6.01652 

0.57227 



11 

1 . 08676 

0.18592 

0.20327 

0.68298 

0.07702 

1.17796 

0.53280 

6.63398 

0.57 376 



12 

1 .06669 

0.20026 

0.19156 

0.68126 

0.08227 

1.17868 

0.53188 

6.88588 

0 . 57557 



13 

1 .06839 

0.21337 

0 . 18096 

0.68016 

0.08791 

1 .17868 

0.53155 

7.36506 

0.57621 



16 

1 .03175 

0.22555 

0.17139 

0 . 68085 

0.09377 

1 .17795 

0.53277 

7.86161 

0 .57381 



15 

1.01662 

0.23705 

0 .16279 

0.68627 

0 .09962 

1 .17591 

0.53620 

8.26717 

0.567 00 



16 

1 .00293 

0.26806 

0.15506 

0.69108 

0.10627 

1 .17226 

0.56231 

8,57999 

0.55682 



1 7 

0.99058 

0.25866 

0 .16806 

0.70167 

0 .10752 

1 .16679 

0.59132 

8.71203 

0.93672 



IS 

0 . 97950 

0.26906 

0 .16176 

0.71586 

0.10811 

1.15966 

0.56296 

8.58797 

0.91312 



1 9 

0 . 96961 

0.27928 

0.13607 

0.73270 

0 .10678 

1.15136 

0.57638 

8.13810 

0 . 68562 



20 

0 . 96085 

0.28968 

0 .13096 

0.75063 

0 .09663 

1.16261 

0.59025 

7.32026 

0.65688 



21 

0 , 95317 

0.29971 

0.12632 

0.76767 

0.08265 

1.13656 

0.60287 

6.13171 

0 .63065 



22 

0 . 96653 

0 .31005 

0.12219 

0.78083 

0.06307 

1.12836 

0.61251 

6.61826 

0.61007 



23 

0 . 96089 

0.32058 

0.11853 

0.78855 

0.03960 

1.12696 

0.61771 

2.87668 

0.39901 



26 

0 .93628 

0 . 33112 

0.11538 

0.78933 

0.01385 

1.12501 

0.61763 

1.00509 

0 . 39918 



25 

0.93266 

0.36262 

0.11257 

0.78239 

-0.00897 

1.12889 

0.61172 

-0,65681 

0.61190 

RADIATING 

LINE 











NUMBER 

5 

1 

1 . 23705 

-0.16958 

0.28106 

0.69605 

0 . 06666 

1.17296 

0.56116 

3.67128 

0.55719 



2 

1 .23228 

-0.12090 

0.28010 

0.69626 

0.06758 

1.17276 

0 . 56167 

3.90938 

0.55656 



3 

1.22100 

-0.07625 

0.27696 

0.69631 

0.05066 

1.17263 

0.56169 

6 .15975 

0 . 55611 




1 . 20657 

-0.03131 

0.26638 

0.69610 

0 .05360 

1.17263 

0.56170 

6.60367 

0 . 55609 



5 

1 . 18668 

0.00710 

0.25565 

0.69568 

0.05665 

1.17271 

0,56155 

6 .65500 

0.55638 




1.16296 

0.06093 

0.26318 

0.69508 

0.05986 

1.17288 

0.56128 

6.92176 

0 . 55692 



7 

1 . 16056 

0.07056 

0.23037 

0.69631 

0.06325 

1.17311 

0.56090 

5.20551 

0.55768 



3 

1 . 11839 

0 .09668 

0.21757 

0.69332 

0 .06685 

1,17362 

0.56036 

5.50717 

0 .55876 



9 

1.09696 

0.11937 

0.20516 

0.69192 

0.07062 

1.17393 

0.53952 

5.82792 

0 . 56062 



10 

1.07660 

0.13976 

0 .19338 

0.69002 

0.07669 

1.17666 

0.53831 

6.16166 

0 . 56286 



11 

1.05765 

0.15816 

0.18236 

0.68780 

0.07839 

1.17553 

0.53686 

6 . 50236 

0.56576 


12 

1.03957 

0.17489 

0.17214 

0.68595 

0.08232 

1.17621 

0.53571 

6.84326 

0.56800 

13 

1.02302 

0.19038 

0.16277 

0.68546 

0.08613 

1.17622 

0.53570 

7.16189 

0.56802 

l<i 

1.00773 

0.20485 

0.15416 

0.68700 

0.08947 

1.17526 

0.53730 

7.42036 

0.56485 

15 

0.99369 

0.21853 

0.14629 

0.69094 

0.09195 

1.17318 

0.54077 

7.58037 

0.55794 

16 

0.98085 

0.23162 

0.13911 

0.69734 

0.09304 

1.16997 

0.54611 

7.59977 

0.54726 

17 

0.96916 

0.24426 

0.13254 

0.70577 

0.09209 

1.16586 

0.55288 

7.43403 

0.53363 

18 

0.95858 

0.25659 

0.12655 

0.71541 

0.08843 

1.16126 

0.56039 

7.04627 

0.51842 

19 

0.94905 

0.26872 

0.12108 

0.72521 

0.08153 

1.15671 

0.56778 

6.41471 

0.50335 

20 

0.94055 

0.28077 

0.11613 

0.73405 

0.07125 

1.15272 

0.57418 

5.54387 

0.49020 

21 

0.93301 

0.29251 

0.11165 

0.74076 

0.05777 

1.14986 

0.57876 

4.45901 

0.48074 

22 

0.92641 

0.30493 

0.10762 

0.74442 

0.04170 

1.14849 

0.58091 

3.20631 

0.47625 

23 

0.92073 

0.31719 

0.10403 

0.74466 

0.02409 

1.14876 

0.58046 

1.85278 

0.47714 

2 ^ 

0.91601 

0.32943 

0.10091 

0.74179 

0.00578 

1.15043 

0.57777 

0.44628 

0.48263 

25 

0.91186 

0.34287 

0.09799 

0.73662 

-0.01340 

1.15307 

0.57356 

-1.04241 

0.49132 

RADIATIHG LINE 










NUMBER 6 










1 

1.16846 

-0.17173 

0.23602 

0.70453 

0.04465 

1.16877 

0.54810 

3.62608 

0.54329 

2 

1.16460 

“0.13088 

0.23528 

0.70441 

0.04816 

1.16872 

0.54818 

3.91097 

0.54311 

3 

1.15641 

-0.09105 

0.23179 

0.70404 

0.05187 

1.16877 

0.54810 

4.21344 

0.54329 

6 

1.14453 

-0.05328 

0.22590 

0.70345 

0.05543 

1.16893 

0.54783 

4.50516 

0.54382 

5 

1.12985 

-0.01823 

0.21815 

0.70259 

0.05892 

1.16921 

0.54737 

4.79367 

0.54475 

6 

1.11331 

0.01383 

0.20913 

0.70138 

0.06241 

1.16966 

0.54662 

5.08508 

0.54625 

7 

1.09573 

0.04292 

0.19936 

0.69972 

0.06588 

1.17033 

0.54552 

5.37895 

0.54847 

8 

1.07778 

0.06926 

0.18930 

0.69748 

0.06930 

1.17127 

0.54396 

5.67431 

0.55158 

9 

1.05996 

0.09316 

0.17930 

0.69468 

0.07260 

1.17248 

0.54194 

5.96613 

0.55561 

10 

1.04261 

0.11493 

0.16957 

0.69172 

0.07575 

1.17376 

0.53980 

6.24916 

0.55987 

11 

1.02595 

0.13491 

0.16028 

0.68920 

0.07875 

1.17483 

0.53802 

6.51892 

0.56341 

12 

1.01014 

0.15340 

0.15151 

0.68784 

0.08151 

1.17533 

0.53718 

6.75773 

0.56509 

13 

0.99525 

0.17066 

0.14330 

0.68807 

0.08371 

1.17510 

0.53758 

6.93662 

0.56429 

14 

0.98134 

0.18692 

0.13564 

0.69002 

0.08515 

1.17406 

0.53931 

7.03454 

0.56085 

15 

0.96841 

0.20240 

0.12852 

0.69350 

0.08543 

1 .17234 

0.54218 

7.02258 

0.55514 

16 

0.95644 

0.21725 

0.12192 

0.69804 

0.08414 

1.17019 

0.54575 

6.87337 

0.54799 

17 

0.94544 

0.23163 

0.11583 

0.70301 

0.08095 

1.16791 

0.54950 

6.56869 

0.54044 

18 

0.93539 

0.24567 

0.11023 

0.70778 

0.07561 

1.16583 

0.55293 

6.09720 

0.53354 

19 

0.92625 

0.25948 

0.10508 

0.71173 

0.06799 

1.16423 

0.55554 

5.45715 

0.52825 

20 

0.91800 

0.27318 

0.10038 

0.71422 

0.05822 

1.16341 

0.55687 

4.65980 

0.52550 

21 

0.91060 

0.28686 

0.09610 

0.71472 

0.04652 

1.16357 

0.55658 

3.72369 

0.52605 

22 

0.90404 

0.30062 

0.09223 

0.71295 

0.03330 

1.16481 

0.55452 

2.67436 

0.53016 

23 

0.89836 

0.31431 

0.08876 

0.70894 

0.01904 

1.16706 

0.55078 

1.53863 

0.53762 

24 

0.89336 

0.32847 

0.08560 

0.70273 

0.00362 

1.17025 

0.54547 

0.29535 

0.54818 

25 

0.88898 

0.34341 

0.08263 

0.69450 

-0.01695 

1.17419 

0.53887 

-1.39848 

0.56129 




OUTPUT 

continues for all 97 radiating lines 



Plot Output 










The graphic output 

(figs 

. 9 to 15) was produced by a local graphics 

package at Lewis 

, by using the plot data 

save file. 






Figure 9. - Solution mesh. 



Figure 10. - Surface Mach number. 
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Figure 15. - Trailing-edge region. 
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APPENDIX - SYMBOLS 

A variable coefficient of change in potential for field equation 

(eqs. (19), (21)); flow area normal to stream-channel direction 
(eqs. (26) to (31)) 

Ag total surface area of finite volume element 

a contravariant velocity component (n direction) in computational 

space 

b contravariant velocity component (c direction) in computational 

space 

C circulation about airfoil as defined by given values of flow angle 

c blade chord 

length of passage along x axis, from leading edge to trailing 
edge when blade is at setting angle 

ex,6y direction cosines, defined in eq. (15) 

G transposed terms in field equation involving known (previously 

determined) potential corrections 

h local stream-tube height normalized to unity far upstream 

i, length of bounding surfaces for finite volume (fig. 4) 

M Mach number relative to blade 

unit vector, normal to airfoil surface in eq. (34); normal to faces 
of finite volume in eqs. (4), (5), and (7) 

q magnitude of flow speed, normalized to unity at far upstream 

conditions 

R residual net mass flow out of a finite volume calculated by using 

potential values from previous iteration (eq. (6) or (16)) 

aR change in net mass flow out of a volume caused by changes 6<p in 

potential values 

^ total residual right side for field equation at each volume, -R + 6 

r local spanwise radial position of stream channel normalized to 

unity far upstream 

S flow direction as defined by contravariant velocities a and b 

in eq. (24) 

s blade spacings 
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X component of relative velocity in physical plane 
change in relative streamwise velocity caused by changes 69 
total absolute velocity vector, uT + vj + sJrj 

volume of element surrounding a grid point; equals area of element 
within stream surface times stream-channel thickness 

y component of relative velocity in physical plane 

change in relative tangential velocity caused by changes 69 

total relative velocity vector, uT + vj 

change in total relative velocity caused by changes 69 

position coordinate along stream channel (fig. 1) 

tangential position coordinate (fig. 1) = ri_E • e 

angle of flow as measured from x axis 

angle of flow as measured from chordline 

ratio of specific heats 

coordinate in computational space whose contours in physical space 
radiate from blade and terminate on periodic boundary 

coordinate in computational space whose physical contours surround 
a i rf 0 i 1 

tangential coordinate in radians, y = 0 • ri^ 

stagger or setting angle of blade row, as measured from x axis to 
chordline 

switch function to control artificial density modification 
(eq. (23)) 

3.14159265 

physical density 

change in physical density caused by changes 69 
artificial density as described in eq. (22) 
flow velocity potential function 

change in local potential field; unknowns in eq. (14) 



n angular velocity of blade row, rad/sec 

0 ) overrelaxation parameter, 1 < w < 2 

Subscripts: 

a absolute frame of reference fixed in space 

av average 

i position on computational plane in ? direction 

j position on computational plane in n direction 

LE at leading edge of blade 

c differentiation with respect to c 

n differentiation with respect to n 

in conditions at far upstream station 

out conditions at far downstream station 

Superscripts: 

k,K artificial time step (iteration level) 

(~) average over the length i of velocity component normal to 

element surface along i (eq. (16) and fig. 4) 
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